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PREDICTING  TRAFFIC  IMPACTS  AT  TWO-LANE  HIGHWAY  WORK  ZONES 


Implementation  Report 

Construction  or  maintenance  activity  on  a  two-lane  highway  often  requires  a  lane 
closure.  As  only  a  single  lane  serves  traffic  in  both  travel  directions,  right-of-way  to 
traverse  the  work  zone  is  sequentially  allocated  to  each  directional  movement  by  flaggers 
stationed  at  both  ends  of  the  work  zone.  Motorist  delays  resulting  from  the 
implementation  of  this  alternating  traffic  control  are  often  significant.  Thus,  a-priori 
evaluation  of  the  resulting  delays  becomes  useful  for  determining  expected  impacts 
created  by  proposed  work  activity  and  for  identifying  appropriate  work  zone  operating 
strategies  such  as  work  zone  physical  length,  hours  of  operation  and/or  the  need  for 
detour  routes  to  minimize  impacts. 

Principles  of  queueing  theory  were  used  to  derive  models  for  estimating  expected 
delay  as  a  function  of  directional  traffic  demand  rates,  work  zone  physical  length  and 
observed  traffic  measures  prevailing  at  work  zones.  The  models  capture  stochastic  effects 
of  expected  delay  as  a  function  of  work  zone  operating  conditions.  To  evaluate  the 
accuracy  of  the  proposed  stochastic  models,  a  computer  simulation  model  was  developed. 
Based  upon  the  outcomes  generated  by  the  simulation  model,  the  delay  models  appear  to 
adequately  predict  the  impacts  of  two-lane  highway  lane  closures. 

To  simplify  and  expedite  the  implementation  of  the  proposed  methodologies,  all 
formulated  procedures  have  been  fully  computerized.  A  user  manual  for  executing  the 
computer  software  is  included  in  the  final  report.  The  software  can  be  run  on  IBM- 
compatible  personal  computers  (versions  286  or  higher)  and  a  math  co-processor  is 
required  on  286  machines.  It  is  a  free  standing  PC  package  using  standard  DOS 
commands. 
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1.  INTRODUCTION 


A  two-lane  highway  is  defined  as  a  two-directional  roadway  with  one  lane  available  for 
traffic  traveling  in  each  direction  [TRB,  1985].  Maintenance,  construction  or  rehabilitation 
activity  on  a  two-lane  highway  often  requires  a  lane  closure.  The  highway's  remaining  lane 
serves  two-directional  traffic  within  the  "work  zone"  by  providing  alternating  right-of-way. 
Flaggers,  stationed  at  both  ends  of  the  work  zone,  sequentially  allocate  right-of-way  to  each 
directional  movement  (vehicles  may  or  may  not  be  constrained  to  follow  a  "pilot  car"  while 
traversing  the  work  zone).  This  operating  scheme,  commonly  referred  to  as  "one-way  traffic 
control",  is  illustrated  in  Figure  1-1. 

Motorist  delays  resulting  from  the  implementation  of  one-way  traffic  control  are  often 
significant  [California,  1978].  Thus,  a  priori  evaluation  of  the  resulting  delays  becomes  useful 
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1)  for  determining  expected  operational  impacts  created  by  proposed  work  activity  and  2)  for 
identifying  appropriate  work  zone  operating  strategies  such  as  work  zone  physical  length,  hours 
of  operation  and/or  the  need  of  detour  routes  to  minimize  impacts. 


1.1.    Previous  Research 

Previous  research  concerning  delay  prediction  under  one-way  traffic  control  has  been 
sparse.  Ceder  and  Regueres  [1990]  used  simulation  to  obtain  the  average  work  zone  delays  and 
compared  these  outcomes  with  the  average  delays  estimated  using  Webster's  delay  model  [1966] 
for  fixed  time  signalized  intersections.  Differences  in  delay  values  yielded  from  the  two  models 
exhibited  a  tendency  to  increase  exponentially  with  decreasing  cycle  time.  To  capture  the 
discrepancy  between  the  two  models,  regression  was  used.  The  regression-based  error  expression 
was  added  to  Webster's  equation.  This  "extended"  model  was  constructed  for  the  special  case  of 
equal  right-of-way  times  for  both  directions  of  travel  and  was  "calibrated"  for  work  zones  of 
relatively  short  physical  length. 

Aside  from  the  obvious  application  limitations  inherent  in  this  approach,  the  fundamental 
soundness  of  the  resulting  prediction  model  is  suspect.  The  regression-based  error  term  developed 
in  this  previous  work  reflects  only  the  characteristics  prevailing  under  the  (narrow)  range  of 
operating  conditions  evaluated  in  the  study.  More  importantly,  Webster's  equation  is  a  steady- 
state  model  for  predicting  delay  at  fixed-time  signalized  intersections,  while  cycle  length  and 
right-of-way  times  are  variable  under  one-way  traffic  control. 

Cassidy  and  Han  [1993]  proposed  a  technique  for  estimating  vehicle  delays  and  queue 
lengths  on  two-lane  highways  operating  under  one-way  traffic  control.  The  model  used 
deterministic  queueing  analysis  techniques  to  formulate  equations  for  computing  average  motorist 
delays  and  queue  lengths  occurring  over  the  given  period.  The  procedure  estimated  average 
required  right-of-way  times  and  relied  upon  empirically  estimated  values  for  traffic  stream 
characteristics  such  as  average  queue  discharge  rate. 
The  model  exploited  a  deterministic  queueing  approach,  which  assumed  uniform  vehicle  arrival 


and  departure  rates  and  used  only  average  values  for  relevant  parameters.  Thus,  the  model  does 
not  consider  the  variability  inherent  in  one-way  traffic  control. 

Cassidy,  Son  and  Rosowsky  [1994]  made  an  effort  to  capture  the  stochastic  nature  of 
work  zone  operations.  The  work  employed  statistical  estimation  and  Monte-Carlo  simulation  to 
identify  distributions  of  average  delay  as  a  function  of  conditions  prevailing  at  the  work  zone. 
Approximate  analysis  techniques  were  then  used  to  compute  user-specified  percentile  values  of 
average  delay  per  operating  sequence  per  direction  of  travel.  Percentile  values  of  average  delay 
were  used  to  probabilistically  evaluate  operating  conditions.  The  approximate  techniques 
produced  easy  to  use,  closed-form  delay  expressions.  However,  the  use  of  Monte-Carlo 
simulation  to  simply  generate  random  variables  has  several  inherent  deficiencies  : 

1.  The  approach  implicitly  assumes  that  a  randomly  arriving  vehicle  is  equally  likely  to 
arrive  in  any  cycle  irrespective  of  cycle  length.  In  fact,  a  randomly  arriving  vehicle  is 
more  likely  to  arrive  during  a  longer  cycle. 

2.  The  variability  in  discharging  headways  exhibited  by  queued  vehicles  will  be  a  function 
of  the  number  in  queue.  This  is  not  reflected  in  the  approach  used  for  generating  values 
of  saturation  flow  rate. 

3.  The  generated  delay  measure  reflect  an  average  value  which  is  exceeded  in  a  specified 
percentage  of  cycles.  This  measure  is  not  equivalent  to  a  percentile  value  of  individual 
vehicle  delay  which  might  better  characterize  work  zone  impacts. 

Most  notably,  the  predicted  delays  reflect  average  values  per  cycle  rather  than  an  expected 
value  over  some  extended  time  period  of  interest.  A  single  expected  value  would  characterize 
operating  conditions  more  definitively.  Moreover,  an  expected  value  would  not  require  the  analyst 
to  subjectively  select  a  priori  some  percentile  value  for  assessing  conditions. 

Although  not  directly  applied  to  two-lane  highway  work  zone  operations,  Newell  [1969] 
used  fluid  and  diffusion  queueing  approximations  to  analyze  the  behavior  of  vehicle-actuated 
signals  at  the  intersection  of  two  one-way  streets  with  no  turning  traffic.  Recognizing  that  one- 


way  traffic  control  operations  behave  in  a  very  similar  manner  to  this  "special  case"  of  a  vehicle- 
actuated  intersection,  Newell' s  procedure  can  be  extended  to  model  two-lane  highway  work 
zones. 


1.2.  Research  Objectives  and  Scope 

The  objective  behind  this  research  has  been  to  formulate  analytical  expressions  (i.e., 
equations)  for  estimating  expected  vehicle  delay  as  a  function  of  work  zone  operating  conditions 
such  as  work  zone  physical  length,  work  activity  type  and  directional  traffic  demands.  The  work 
has  employed  fluid  queueing  approximations  previously  developed  by  Newell  [1969]  for 
estimating  delay  at  vehicle-actuated  intersections  of  two  one-way  streets.  The  model  captures 
stochastic  affects  on  expected  delay  as  a  function  of  directional  traffic  demands,  work  zone 
physical  length  and  observed  traffic  stream  measures  prevailing  at  under-saturated  work  zone 
operations.  The  term  under-saturated  denotes  that  directional  queues  are  completely  served  during 
each  allocation  of  right-of-way. 

For  over-saturated  conditions  (i.e.,  operations  where  residual  queueing  prevails), 
deterministic  queueing  techniques  are  used  for  predicting  average  delays  over  a  specified  time 
interval  of  interest  [Cassidy,  Son  and  Rosowsky,  1993]. 

To  render  the  proposed  techniques  simple  and  straightforward  to  use,  the  models  have 
been  computerized  as  part  of  a  user-friendly  software  package.  As  such,  application  of  the 
derived  procedures  is  relatively  simple  and  timely.  A  "User  Manual"  for  this  software  is 
contained  in  Appendix  B  of  this  report. 

To  insure  that  the  resulting  predictions  accurately  estimate  delay,  model  predictions  have 
been  carefully  validated  using  both  empirical  measurements  and  the  outcomes  generated  from  a 
microscopic  simulation  model  developed  as  part  of  this  research  project. 


To  insure  broad  applicability,  delay  prediction  techniques  have  been  developed  for  a 
number  of  differing  work  activity  types.  These  activity  types  are  : 

•  asphalt-concrete  overlays  (with  pilot  car  operation) 

•  chip-seal  resurfacing  (with  pilot  car  operation) 

•  general  construction  (i.e.,  shoulder  work  ,  saw  cutting,  pavement  marking)  without  pilot 

car 

•  general  construction  with  pilot  car 

Where  a  work  activity  type  not  addressed  above  is  to  be  analyzed,  approximate  predictions  can 
be  generated  by  specifying  one  of  the  four  work  activity  types  above  that  most  closely  "matches" 
the  activity  in  question. 


1.3.  Report  Overview 

The  models  derived  for  each  activity  type  are  generally  identical.  Major  differences 
between  the  models  are  reflected  only  in  the  distributions  of  prevailing  traffic  stream  parameters. 
For  example,  the  queue  discharge  rates  and  work  zone  travel  speeds  under  asphalt-concrete 
overlay  (A-C)  activities  are  significantly  different  from  those  prevailing  under  chip-seal 
operations. 

This  report  therefore  details  all  tasks  involved  in  developing  models  for  evaluating 
asphalt-concrete  overlay  operations.  The  relevant  operational  parameters  under  all  other  work 
activity  types  are  summarized  in  Appendix  A. 

Chapter  2  of  this  report  presents  the  derivations  of  the  stochastic  queueing  model  to 
estimate  work  zone  operating  characteristics  and  the  resulting  directional  average  delays.  Included 
in  chapter  2  is  a  description  of  the  queueing-based  assessment  technique  for  over-saturated 
conditions.  Chapter  3  describes  empirical  data  collection  and  analysis.  Chapter  4  describes 


development  of  the  microscopic  simulation  model  and  the  simulation  experiments  associated  with 
validating  the  proposed  procedures. 

A  summary  and  conclusions  are  presented  in  the  fifth  chapter  of  this  report.  Appendix  A 
presents  parameter  values  for  work  activity  other  than  A-C  operations.  Appendix  B  presents  a 
"User  Manual"  for  executing  the  computer  software  developed  in  this  project. 


2.  QUEUEING  DELAY  MODEL 


2.1.    Introduction 

Given  that  the  implementation  of  one-way  traffic  control  typically  results  in  significant 
motorist  delay  [California,  1978],  a-priori  evaluation  of  work  zone  operation  should  be  performed 
to  identify  expected  impacts  and/or  to  establish  suitable  work  zone  strategies  (i.e.,  the  physical 
length  of  the  lane  closure,  hours  of  work  zone  operation,  the  need  for  diversion  routes,  etc). 

The  work  presented  in  this  chapter  has  used  principles  of  queueing  theory  to  derive 
closed-form  expressions  for  estimating  expected  delay  as  a  function  of  work  zone  operating 
conditions.  The  work  has  exploited  fluid  approximations  and  stochastic  models  originally 
developed  by  G.F.  Newell  [1969]  for  predicting  delay  at  the  intersection  of  two  one-way  streets 
controlled  by  a  vehicle-actuated  traffic  signal.  Traffic  stream  characteristics  under  one-way  traffic 
control  are  similar  to  those  prevailing  at  the  stop  bar  of  a  signalized  intersection  [Cassidy  and 
Han,  1993;  Ceder  and  Rogueros  1990].  Thus,  signalized  intersection  delay  models  serve  as  a 
logical  foundation  for  work  zone  delay  prediction  models. 

The  primary  difference  between  traffic  signal  and  work  zone  operation  is  the  significance 
of  lost  time.  For  traffic  signals,  clearance  intervals  (i.e.,  yellow  times)  are  fixed  and  account  for 
a  small  portion  of  the  overall  cycle  length.  In  contrast,  the  lost  times  occurring  under  one-way 
traffic  control  include  platoon  travel  times  through  extended  work  zone  sections  as  well  as  the 
lost  times  prevailing  at  either  end  of  the  lane  closure  during  the  initiation  of  directional  right-of- 
way.  As  lost  time  accounts  for  a  significant  portion  of  the  cycle  length  and  is  stochastic  in 
nature,  Newell' s  delay  expressions  for  signalized  intersections  were  altered  to  capture  these 
effects. 

Work  zone  operation  also  differs  from  typical  signalized  intersection  operation  in  that 
flaggers  commonly  extend  right-of-way  well  beyond  queue  dissipation  to  accommodate  arriving 
vehicles.  Thus,  expressions  were  derived  to  account  for  the  influence  of  large  green  time 
extensions   (although   inclusion   of  these   effects   were   ultimately   found   to   provide   little 


improvement  in  model  accuracy). 

Recognizing  again  that  queueing  characteristics  at  a  work  zone  entrance  are  similar  to 
those  at  the  stop  bar  of  a  signalized  intersection,  terminology  typically  associated  with  signal 
operation  have  been  adopted  for  describing  traffic  flow  characteristics  under  one-way  traffic 
control. 

•  Effective  Green  Time,  G,  is  the  right-of-way  time  utilized  by  a  single  directional 
movement  in  one  cycle.  It  is  the  elapsed  time  between  the  first  vehicle  and  last  vehicle 
entry,  traveling  in  the  same  direction,  during  a  single  allocation  of  right-of  way. 

•  Cycle  Length,  C,  is  the  time  required  for  one  complete  operational  sequence.  For 
example,  cycle  length  may  be  defined  as  the  elapsed  time  between  successive  allocations 
of  effective  green  time  for  a  single  directional  movement. 

•  Saturation  Green  Time,  g,  reflects  the  portion  of  the  effective  green  in  which  queued 
vehicles  are  discharged  into  the  work  zone.  It  is  the  elapsed  time  between  the  first 
vehicle  and  last  queued  vehicle  entry. 

•  Saturation  Flow  Rate,  S,  reflects  queue  discharge  rate  into  the  work  zone  during  saturated 

green  time. 

•  Start-Up  Lost  Time,  L,  is  defined  as  the  elapsed  time  between  the  last  vehicle  in  the 
opposing  direction  exiting  the  work  zone  and  the  entry  of  the  first  queued  vehicle 
traveling  in  the  subject  direction. 

•  Green  Time  Extension,  E,  is  the  unsaturated  green  time.  It  reflects  that  portion  of  the 
effective  green  time  extending  from  queue  dissipation  to  the  final  vehicle  entry  into  the 
work  zone.  Discharge  rate  during  this  period  is  equal  to  vehicle  arrival  rate. 


Effective  Red  Time,  R,  is  the  amount  of  time  per  cycle  that  effective  green  time  is  not 
allocated  to  a  given  directional  movement. 


2.2.   Fluid  Approximations 

To  illustrate  the  characteristics  of  one-way  traffic  control,  we  begin  by  evaluating 
operation  using  fluid  approximations.  The  evolution  of  opposing  vehicle  queues  at  the  work  zone 
is  depicted  in  Figure  2-l(a).  During  cycle  j,  traffic,  presumed  to  behave  like  a  fluid,  arrives  in 
direction  i  at  rate  qt,  where  i  =  1  or  2.  A  queue  in  direction  i  builds  during  the  red  time  of  the 
y'th  cycle,  Rir  During  the  "saturated"  green  time  of  the  y'th  cycle  in  direction  i,  gt],  vehicles 
discharge  into  the  work  zone  at  rate  5,. 

Figure  2- 1(a)  is  identical  to  the  queueing  diagram  for  the  intersection  of  two  one-way 
streets  controlled  by  a  vehicle-actuated  (V-A)  signal  presented  in  Figure  2- 1(b).  The  two  figures 
exhibit  only  minor  differences.  For  work  zone  operation,  values  of  Rtj  are  not  only  functions  of 
the  saturated  green  time  in  the  opposing  direction,  but  also  the  travel  time  of  both  directional 
platoons  through  the  work  zone  and  an  added  lost  time  at  the  opposing  end  of  the  work  zone. 
We  define  Ttj  as  the  amount  of  time  for  the  rear  of  a  platoon  to  traverse  the  work  zone  in 
direction  i  during  the  y'th  cycle.  Values  of  T*  may  well  vary  by  direction  [Ceder,1993]  while 
values  of  L  are  small  relative  to  cycle  length  and  are  therefore  treated  deterministically  and 
independent  of  travel  direction. 

Finally,  directional  effective  green  time,  G{j,  is  often  extended  by  flaggers  at  the  work 
zone  entrance  well  beyond  queue  dissipation  to  accommodate  arriving  vehicles.   Thus, 

G    =  g..  +  E  (2.1) 

where: 

Etj        =  green  time  extension  beyond  queue  dissipation  in  direction  i  during  y'th 

cycle 
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Figure  2-1  (a) 

Queue  Length  vs  Time 
(Two-Lane  Highway  Work  Zone) 
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Figure  2-  1(b) 

Queue  Length  vs  Time 
(Vehicle-Actuated  Intersection  -  One  Way  Street) 
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Vehicles  in  each  direction  are  treated  like  a  fluid  with  a  uniform  arrival  rate  of  g,  and 
a  queue  discharge  rate  of  5,. 

Defining  t=0  as  the  end  of  G^.,  (i.e.,  the  end  of  effective  green  time  in  direction  i  during 
the  j-hh  cycle)  and  Qt(t)  as  the  queue  length  at  time  t  in  direction  i, 

Off)  =qt-t  o  <  t  <  t R 

=  q.-t  -  5,  •  (t-tR)  tR<t<tg  (2.2) 

=  0  t   <  t  <  tr 


where  tR,  tg  and  tG  represent  the  end  times  of  the  red,  saturated  green  and  total  green  periods, 
respectively. 

Assuming  Ri}  is  known  and  noting  that  directional  queues  do  not  increase  during  their 
respective  green  periods, 

«i  <?V  +  82j  -  E2j  *  T2j  +  2L)  =  (Srqi)glj  (2.3) 

«2  CTV  -  gy  -  EXj  +  T1M  +  2L)  =  (S2-q2)g2M  (2.4) 

By  defining  lost  time,  Yip  as  the  portion  of  cycle  j  where  saturated  green  periods  are  not 
occurring  for  either  travel  direction, 


Y,.  =  rr  +  £,.  +  r,.  +  2L  (2.5) 

U  U  ?/  2/ 


K.  =  r,.  +  £..  +  r...  +  2L  (2.6) 

2/  2/  I7  ly+1 


Directional  red  time,  Rip  is  the  sum  of  total  directional  lost  time  and  the  saturated  green 
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period  for  the  opposing  direction. 

*„  =  S2J  +  Yv  (2.7) 

**,  =  gv  +  Y2J  (2-8) 

From  Equations  (2.3)  and  (2.4),  the  saturated  green  times  become 


*v  -  7J    f  (2-9) 


(5,  - <?,) 

=     ^Ij     +    Yl)  (2   1Q) 

(5,   -  q2) 


For  Equations  (2.9)  and  (2.10),  the  saturated  directional  green  time  during  a  given  cycle 
can  be  evaluated  from  its  value  in  the  previous  cycle,  i.e., 

«A  (Y-Y)+^Y 


8ij+\ 


qjlgy  +    S,S2       Xl     y       S2     *  (2.11) 

(5,-g,)  (52-<?2)   +     (1  -  g/S^Cl  -  <?A) 


Assuming  work  zone  operation  is  under-saturated  (i.e.,  residual  queueing  does  not  occur), 
the  fluid  approximation  infers  that  asy— »  °°,  Yy— >  y,   ,Etj  — >  Ei  and  gtj  — >  g,  [Newell,  1969]. 

5,S2  Ul     2)       S2     2  (2.12) 

*2       (i  -  <? A  -  9 A) 


The  subscripts  in  Equation  (2.12)  can  be  reversed  for  estimating  g;.  If  Y,  =  Y2,  Equation  (2.12) 
is  the  same  as  the  formula  derived  by  Newell  and    g2j  can  be  expressed  as 
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Sn  =  S2  +  (  g2l  -  82  )  «;M 


>2j  °2 


(2.13) 


Where 


a  = 


q& 


\ii 


{sx-qx)(s2-q2) 


Finally,  cycle  length,  C,  becomes  the  sum  of  Y,,  g;  ,  g2  and  E, 


C*Y.+E.+g.+* 


=  Y    +  E    + 

1  j  A^j 


5ly  +  ±y 

(1  -  <?  A  -  q A) 


(2.14) 


The  total  delay  for  all  queued  vehicles  in  direction  1  during  time  interval  t  to  t+dt  is 
Q,(t)dt.   Thus,  the  total  directional  delay  per  cycle,  Wu  becomes 


w,  =  ry,%  ceo* 


<?,  (fe  +  rt)2 

2   (i  -  <?A>) 


(2.15) 


Once  again,  delays  for  the  opposing  direction  can  be  estimated  by  reversing  the  subscripts. 


And  the  average  delay  per  vehicle  for  both  directions,  w.  ,  is 


<7,  q 


1    12 


W.    = 


s,s2 


4\    4i. 


2(q+q2){\-lL-lL)  \K 


( 


{yx-y2Y-Lyx 

12 


f 


V<7i) 


\ 


(y2-r,)+_Ly2  [s2-q2) 


(2.16) 


sx     S2 
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If  one  were  to  assume  Y,  =  Y2,  an  inappropriate  presumption  for  work  zones  according 
to  Ceder  [1993],  the  resulting  expression  for  average  vehicle  delay  in  both  directions  becomes 
identical  to  the  expression  derived  by  Newell  [1969]  for  V-A  intersections. 


2.3.   Stochastic  Modeling 

The  delay  expression  reflecting  fluid  approximations  (Equation  (2.15))  does  not  account 
for  fluctuations  in  vehicle  arrival  and  departure  processes  and  other  stochastic  properties  of  work 
zone  operation.  Figure  2-2  depicts  the  possible  realizations  of  directional  queue  length  versus 
time,  Q,{t)  [Darroch,  Newell  and  Morris,  1964].  The  solid  line  in  Figure  2-2  represents  the 
expected  evolution  of  queue  length  conditioned  on  a  known  red  time.  As  expressed  in  Equation 
(2.15), 


E(W\R.)  =  il        (^)2  (2.17) 

''    ljJ        2    (1   -  q/S) 


where  the  subscript  j  is  used  in  Equation  (2.17)  to  denote  the  cyclic  variability  of  parameters. 

Treating  directional  vehicle  arrivals  to  time  t,  At(t),  as  a  stationary  process 

E  \A0  =  q:  t  (2.18) 

We  exploit  Newell' s  assumption  that  there  is  a  variance  to  mean  ratio  of  the  arrival  process,  Iia, 
of  the  order  1  that  is  independent  of  the  mean. 

War  IA0  =  Laqt  (2.19) 

We  also  define  the  departure  process  during  the  saturated  green  time,  £>/?),  identically  to  Newell. 
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Q(t)IR 


\Q'(t)IR 
J(Q(t)IR^\\ 

\  \  \ 


E(glR) 


Figure  2-2 
Evolution  of  Random  Queue 


Namely, 


E  <D0  =  Sfr-a) 


(2.20) 


where  a,  is  a  correction  due  to  "start-up"  delays  immediately  following  the  initiation  of  the  green 
period  and  t=0  at  this  green  initiation.  The  variance  of  directional  vehicle  departure  time  t  is 


Var  \D0  =  Vi  +  IJ5& 


(2.21) 


where  the  variance  to  mean  ratio  Iid  is  again  of  order  1 .   Start-up  delays  can  be  captured  in  the 
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term  V,  .  Given  the  very  long  cycle  times  expected,  Vt  is  small  relative  to  the  second  term  in 
Equation  (2.21). 

The  dashed  lines  in  Figure  2-2  denote  two  possible  realizations  of  the  random  queue 

length  Qt  (0  (one  above  the  average  and  one  below)  caused  by  fluctuations  in  vehicle  arrival 

and/or  departure  rates.  If  the  queue  dissipates  before  the  end  of  the  expected  saturated  green  time, 
the  queue  could  be  negative  at  the  end  of  this  time.  If  the  queue  does  not  dissipate,  the  queue 
is  positive  at  the  end  of  expected  saturated  green  time.  Whether  the  directional  queue  in  cycle 
j  dissipates  before  or  after  the  end  of  the  expected  saturated  green  period,  the  shaded  triangular 
area  representing  the  added  delay  (which  may  be  positive  or  negative)  in  Figure  2-2  can  be 
approximated  knowing  that  the  slope  of  Qt(t)  during  the  saturated  green  time  is  approximately 
(qrSj).  Thus,  the  "extra"  directional  delay  caused  by  fluctuating  arrival  and  departure  processes, 
Wie,  becomes 

W.lR-l   gjgW  +  *»»  (2-22) 

bl    ij       2  (S.  -  q) 

where  queue  length  at  the  end  of  the  expected  saturated  green  period  is  the  difference  between 
cumulative  arrivals  and  departures, 

<2,  (  fiCgJjy  +  R,  )  -  A,  (  EQg/y  +  Rt)  )-Di(  Efejfy  *  «,)  )         (2-23) 

Taking  the  expectation  of  Equation  (2.22), 

E(W  \R  )  =  —  J     l]  ,J  (2.24) 

2  (S.  -  q) 

The  expected  number  of  vehicles  in  queue  at  the  end  of  the  expected  saturation  green  period  is 
0,  thus 
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E(  QiXE(g..\Ry+R.p)  ) 


Var(Q.(E(g.j+Rip)) 


=  VariAfMg,  \Rt)  +Rtj))  *V*Q>f&Mi  Ity))  (2.25) 

IqR 

=       '  '   »      +  V 


Adding  Equation  (2.24)  to  the  expected  delay  as  expressed  by  the  fluid  approximation 
(Equation  (2.17)),  the  total  expected  delay  due  to  variation  in  the  arrival  and  departure  processes 
becomes 


q.  R- 

E(W..\R.)  =         '     ' 
"'    +       2(1 -q/S) 


1   + 


/ 


V. 


S{l-qWRv       Sfltf 


(2.26) 


where    /,  =  lia  +  Iid 


Taking  the  expectation  over  the  distribution  of  Rip 


E(W.)  =  - 

y        2(1 -q/S) 


E\Rt)  +  VariRJ 


S^l-q/S)        Sfli 


(2.27) 


Equations  (2.26)  and  (2.27)  are  identical  to  expressions  previously  derived  for  signalized 
intersections  [Newell,  1969].  At  this  point,  we  adapt  the  intersection  delay  equations  above  to 
account  for  the  operating  characteristics  unique  to  one-way  traffic  control  at  two-lane  highway 
work  zones.  Namely,  the  stochastic  affects  on  delay  extend  beyond  fluctuating  arrivals  and 
departures.  Variability  in  the  times  required  to  process  platoons  through  the  work  zone,  Ttj,  and 
the  correlated  influence  of  green  time  extensions,  E«,  influence  vehicle  delays  (although  the  latter 
influence  was  ultimately  found  to  be  insignificant). 

The  expected  red  time  in  direction  1  during  the  /th  cycle  is  the  sum  of  several  previously 
defined  random  variables. 
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E(RJ  =    E{  Tv  +  g2.  +  £^  +  r2.  +  2L  ) 

(2.28) 

=  E{TX)  +  £(g2;  +  £(£2;  +  £(r2;  +  2L 

Given  that  each  of  the  n  random  variables  composing  Ru ,  Xt,  have  a  finite  second  moment, 

Var(i  X)  =  Z   Var  X   +  22  2    Cov(X.  ,  X)  (2-29) 

Evaluation  of  a  large  empirical  data  base  [Cassidy,  Son  and  Rosowsky,  1993]  indicated 
that,  with  the  exception  of  T1}  and  g2j,  the  random  variables  in  Equation  (2.28)  are  generally 
uncorrelated  or  only  slightly  correlated.  Added  time  to  "process"  platoons  through  the  work  zone, 
Tjj  does  increase  queue  length  in  the  opposing  direction  (and,  thus,  g2]).  Therefore, 


Cov(g2j  ,  7y  =  E(Tt/  g2j)  -  £(Ty)  •  E(g2j) 


(2.30) 


S2~92 


Var(Tu) 


w 


Given  that  lost  time,  L,  is  treated  deterministically, 


VartRJ  =  if!!*-  VariTJ  *  Var(g2j)  *  Var(E2j)  +  Var(T2)  (2.31) 


2.4.   Estimation  of  Moments 

To  exploit  the  expression  for  expected  directional  delay  (Equation  (2.27)),  the  first  and 
second  moments  of  gtj  and  Etj  are  required. 
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2. 4.1. Saturated  Green  Time 

By  exploiting  previously  used  assumptions  concerning  the  independence  of  arrival  and 
departure  processes  from  one  cycle  to  the  next  [Darroch,  Newell  and  Morris,  1964;  Dunne  and 
Potts,  1965;  Newell,  1969],  the  conditional  mean  saturated  green  times  given  Rtj  derived  by  the 
fluid  approximation  are  approximately  equal  to  gir   That  is,  E(g.p  =  £(£(g  |J?p) .   Thus, 


E(§1J)  =  -IL.  [E(Th)  *  E{g2)  *  £(£2;)  +  E(T2p  +  2L] 


(2.32) 


E(82M)  =  -r^-  [^  +  %v)  +  mx)  +  E(Tlj+l)  +  2L] 

d2  <?2 


(2.33) 


With  a  stationary  process, 


£(gl)  =  J!i_  [ECTJ  +  £(g2)  -  £(£2)  +  £(r2)  +  2L] 

5r4, 


(2.34) 


E(g2)  =  _!L_  [E(T2)  +  £(gl)  +  £(£,)  +  £(7,)  +  2L] 


(2.35) 


From  Equations  (2.34)  and  (2.35),  the  expected  saturated  green  time  (for  direction  1)  is 


E(gl)  = 


qJS 


V"i 


^-^Srq2/S2 


EiTJ+EiTJ+SiEiEjfi-!^.)  £(£2)+2L 


(2.36) 
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where  subscripts  are  reversed  to  estimate  E(g2). 

As  the  saturated  green  time  given  the  red  duration  is  the  sum  of  the  conditional  expected 
value  of  gy  and  the  additional  term  accounting  for  fluctuations  in  arrivals  and  departures  [Newell, 
1969], 


[  A.(  R    +  E(g..\R.)  -  D(  E(R.)  )  ] 

g..\R..  -  E(g..\R.)  +         '      ,J ';     ,J '  ,; 

y     »  y     ,}  (S.-q) 


(2.37) 


The  conditional  variance  of  saturated  green  time  can  be  expressed  as 


v*ig/RJ 


=  Var 


[  AJL  Ry  *  E(gy\Ry)  )  -  D.(  E(R(p  )  ] 


(S-q) 


InR 

1*1       IJ 


5. 


(V<7,)2 


(2.38) 


and  the  unconditional  variance  becomes 


Var(gip  =  E{VaHgJR.))  *  Var(E(g./R.)) 


(2.39) 


Thus,  for  one-way  traffic  control  operations,  again  assuming  a  stationary  process 
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Var(g.) 


S?  (1  -  qJStf 


V, 


(Srqx)2 


(Sx-qx)2 


Var(T{) 


S2-q2 


+  Var(g,)+Var(E2)+Var(T2) 


(2.40) 


Var(g2)  = 


I2q2E(T2  +  gl+El+Tl+  2L) 


V, 


5/  (1  -  qJS2) 


(S2~q2) 


q2 


(S2~q2)2 


Var(T2) 


Var(«i)+Var(£1)+Var(7',) 


(2.41) 


Combining  Equations  (2.40)  and  (2.41),  the  variance  of  saturated  green  time  for  a  given 
direction  is  expressed  without  reference  to  the  saturated  green  time  for  the  opposing  direction. 


2.4.2.  Green  Time  Extension 

We  assume  flaggers  extend  directional  right-of-way  beyond  queue  dissipation  whenever 
the  arriving  vehicle's  headway  is  not  greater  than  some  value  H  (where  the  value  of  H  can  be 
estimated  empirically  using  discrete  choice  techniques  described  in  section  4.1.2).  Hence,  the 
decision  to  extend  directional  green  time  is  a  function  of  the  vehicle  arrival  distribution. 

As  vehicle  inter-arrival  times  can  be  assumed  exponential,  independent,  identically 
distributed  random  variables, 
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E(Xk)  =  1  (2.42) 


where: 


Xk         =  k*  inter-arrival  time 

a,  =  directional  flow  rate 


Letting  the  random  variable  N  be  the  number  of  vehicles  discharged  during  the  green  time 
extension,  and  given  N=n,  the  length  of  the  extended  green  time  is  X,  +  X2  +  ...  +  Xn  and 
E[E.\N=n]  =  E(Xl  +  X2  +  ...  +  X)  =  n  E{XX)  =  nlqr  For  each  event  (i.e.,  vehicle  arrival),  the 

probability  that  the  headway  is  less  than  or  equal  to  H  is  P(  H  >  xk  )  =  e'qH .  Random  variable 

N  has  a  geometric  distribution  with  parameter  p  =  e  ~q'     and  an  expected  value  of  (l-p)/p. 


E(E)  =  E 


N-l 


l-g'g,g.l  (2.43) 


e^     q, 


The  variance  of  Ei  is  also  obtained  by  conditioning  on  the  number  of  vehicles  discharged 
during  the  green  time  extension. 


Var(E)  =  EiVaiiE^N))  +  Var(£(£.|A0) 


\-e-qH    1         \-e-qfI    1 

• +  - 


(2.44) 


-qH  2  -2qH  2 
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2.4.3.  Travel  Time 


Values  of  work  zone  travel  time,  Tt,  clearly  vary  as  a  function  of  work  zone  physical 
length,  C.  Thus,  the  speed  of  the  final  vehicle  through  a  work  zone,  v.,  becomes  the 
generalizable  measure. 


(     \  (     \ 


E(T)  =  E 


v 

V    V 


=  l-E 


1 


V 

^  'J 


(2.45) 


The  moments  of  vr.,  which  vary  as  a  function  of  work  activity  type  [Cassidy,  Son  and 
Rosowsky,  1993],  can  be  empirically  measured. 

Expanding  the  term  —  in  a  Taylor  series   with  constant  C, 
v 


A0  +  A^v.-Q  +  A2(vrC)2  +  Ap-Q'  * 


(2.46) 


A  is  evaluated  by  differentiating  with  j  =  1,2,3,- 


a  -  UL  (1)  ic 

1       7'!  dXJ     V. 


(2.47) 


Letting  C  =  E(vJ  and  taking  the  expectation, 


E(±)=.    '  » 


V.       E(V)     £2(v.) 


E(v-E(v))+^—E(vrE(v))2--J—E(v-E(v))3+"   (2.48) 


£3(v,.) 


E\v) 


By  definition,  £,(v_  -  £(v(.))=0  .  Similarly,  for  any  distribution  that  is  approximately  symmetric, 
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E(vt  -  E(v))a  «  0  with  odd  positive  integer  a.  Thus, 


pfl      E2(v^Variyt) 
V.  E\v) 


(2.49) 


And  the  expected  directional  travel  time  of  the  final  vehicle  through  the  work  zone  becomes, 

£2(v,)  +  Var(y) 


E(T)~  I 


E\v) 


(2.50) 


The  variance  of  travel  time  is  expressed  as 


Var(T) 


=  f  •  Var(-) 
V. 

( 


-    1)2 


=  e 


\ 


i  -,    ->  i 

£(_L)2-£2(_L) 
v  v 


v 


J 


(2.51) 


Expanding  the  term 


(     Y 

1 

v 


in  a  Taylor's  series  and  taking  the  expectation,  we  determine 


i  £2(v.)  +  3Var(v.) 

£(_)2  -  ' _Z 

V,  £*(V.) 


(2.52) 


i  E2(v)  +  3Var(v.) 

Var  (  —  )  «  

'  v.  £4(v,.) 


^^vp+VaKv^ 


£3(v,) 


(2.53) 
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2.5.  Extending  Deterministic  Models  for  Oversaturated  Condition 

Applying  the  above  models  to  evaluate  oversaturated  conditions  at  a  work  zone  entrance 
is  relatively  simple.  The  basic  phenomena  can  be  represented  by  a  queueing  diagram  illustrated 
in  Figure  2-3. 

For  a  known  or  assumed  directional  arrival  flow,  the  cumulative  arrivals  to  time  t,  A(t), 
can  be  constructed  as  shown  in  Figure  2-3.  The  dashed  line  in  Figure  2-3  represents  the  capacity 
of  the  subject  work  zone  entrance.  Its  slope  is  equal  to  (G/C)s  [Hurdle,  1988].  This  dashed  line 
has  been  drawn  tangent  to  Aft)  at  time  zero,  signifying  the  time  when  the  work  zone  entrance 
first  becomes  oversaturated.  The  departure  curve,  Dft),  in  Figure  2-3  represents  the  time- 
dependent  rate  at  which  vehicles  enter  the  work  zone.  Of  primary  importance  is  the  recognition 
that  the  area  between  the  curves  Aft)  and  Dft)  represents  the  total  vehicle  delay  at  the  work  zone 
entrance.  Average  delay  is  this  total  value  divided  by  Aft)  -  AfO).  The  dashed  line  in  Figure  2-3 
thus  divides  the  area  (i.e.,  the  total  delay)  into  two  components  :  uniform  delay  below  the  dashed 
line  and  overflow  delay  above. 

Equation  (2-12)  is  used  to  identify  if  and  when  oversaturated  conditions  prevail  at  a  work 
zone  entrance.  Where  (2-12)  yields  a  required  green  time,  G„  larger  than  permissible  by 
policy/judgement,  the  analyst  is  then  required  to  specify  the  actual  right-of-way  time  which  will 
prevail  at  the  work  zone  entrance.  This  allowable  green  time,  Ga,  should  reflect  the  largest  value 
judged  to  be  acceptable  (i.e.,  the  largest  time  interval  over  which  flaggers  are  willing  to  extend 
right-of-way).  The  value  of  Ga  will  dictate  the  horizontal  distance  (i.e.,  the  time  interval  on  the 
queueing  diagram)  over  which  queue  discharge  will  occur. 

The  effective  red  time  (i.e.,  the  horizontal  distance  over  which  vehicles  do  not  enter  the 
work  zone)  is  the  sum  of  the  required/actual  green  time  in  the  opposing  direction,  the  start-up 
lost  times  and  one  "round-trip"  travel  time  through  the  work  zone. 

To  include  the  delay  component  created  by  reduced  vehicle  speeds  within  the  work  zone, 
the  cumulative  number  of  departing  vehicles  within  the  time  interval  of  interest  is  observed  from 
the  queueing  diagram.  The  product  of  this  value  of  Dft)  and  the  average  travel  time  delay  (i.e., 
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the  difference  between  "actual"  and  "normal"  work  zone  travel  times)  is  added  to  the  total 
queueing  delay. 

The  fact  that  this  analytical  approach  does  not  account  for  delays  created  by  "random 
effects"  will  generally  not  be  a  concern.  Added  delays- created  by  operational  stochasticity  will 
typically  be  very  small  relative  to  the  deterministic  delay  estimate. 

To  eliminate  the  tedium  created  by  requiring  the  user  to  evaluate  delay  by  constructing 
an  appropriate  queueing  diagram,  the  methodology  has  been  computerized  (Appendix  B). 


lime,  t 


Figure  2-3 

Queueing  Diagram  Depicting  Over-saturated  Condition 

(Source  :  Hurdle,  1988) 
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3.  EMPIRICAL  OBSERVATIONS 


The  delay  prediction  models  described  in  the  previous  chapter  rely  upon  empirical 
estimates  of  arrival  rates,  work  zone  travel  speeds,  queue  discharge  rates,  start-up  lost  time  and 
green  time  extensions.  The  variables  qt  and  q0  representing  work  zone  traffic  demands  in  each 
direction  and  work  zone  length  ft  must  be  known  or  assumed  input  to  the  prediction  model.  To 
insure  that  derived  methodologies  accommodate  a  broad  range  of  application,  empirical  data  were 
collected  from  a  number  of  different  work  zone  locations  operating  under  a  variety  of  work  zone 
activity  types.  The  work  has  exploited  a  large  empirical  data  base  collected  in  California  as  part 
of  an  earlier  research  project  sponsored  by  the  California  Department  of  Transportation. 

To  some  extent,  traffic  operations  at  two-lane  highway  work  zones  represent  idealized 
scenarios.  For  example,  traffic  moves  through  a  single  channel  (i.e.,  lane).  Thus,  virtually  no 
vehicle  over-taking  maneuvers  occur.  Moreover,  the  work  zone  areas  used  for  data  collection  in 
this  research  typically  had  few,  if  any,  and  minor  access  and  egress  points.  Thus,  for  all  practical 
purposes,  conservation  of  vehicles  was  maintained  (i.e.,  vehicles  entering  and  exiting  the  work 
zone  were  essentially  the  same).  The  data  collection  techniques  employed  in  the  research  have 
sought  to  capitalize  on  these  idealized  prevailing  conditions.  Chapter  3  summarizes  data 
collection  and  analysis  tasks  to  derive  suitable  estimates.  Table  3-1  summarizes  the  sites  from 
which  data  were  collected.  In  total,  over  70  hours  of  operational  data  were  collected  from 
California  test  sites  reflecting  4  types  of  work  activity. 


3.1.  Data  Collection 

The  data  collection  team  consisted  of  two  persons,  each  equipped  with  a  lap-top  computer 
with  synchronized  internal  clocks.  A  program  (written  in  BASIC  language)  was  developed  to 

record  the  actual  times  in  which  computer  keys  were  pressed.  Specific  keys  of  the  lap-top 
computer  were  designated  for  recording  the  occurrence  times  of  specific  events.  The  resulting 
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Table  3-1 
Summary  of  Empirical  Observation 


SITS 
No. 

LOCATION 

OBSERVATION 
DATS(S) 

WORK 
ACTIVITY 

WORK  ZONE 
LENGTH ( S ) 

PILOT 
CAR 

HOURS  of: 
OBSERV. 

1* 

Highway  28 
Carnelian  Bay 

6/17/91 

3,031  ft 

Yes 

6 

2* 

Highway  23 
Tahoe  City 

6/19/91 

A-C 

5,070  ft 

Yes 

3 

3 

Highway  89 
Truckee 

6/20/91 

A-C 

4,181  ft 

Yes 

3.5 

4 

Highway  193 

Newcastle 

6/21/91 

Pavement 
Markers 

5,964  ft 

No 

0.5 

5 

Highway  267 
Placer  Cnty 

6/24/91 

A-C 

6,481  ft 

Yes 

7 

6 

Highway  193 
Georgetown 

7/8/91  - 
7/10/91 

Chip -Seal 

1,49  5  ft 
-  4.15  mi 

Yes 

11.5 

7 

Highway  267 
Truckee 

7/10/91 

A-C 

5,691 

Yes 

1.5 

8 

Highway  49 

7/15/91  - 
7/16/91 

Chip-Seal 

3.1  - 
3  .4  mi 

Yes 

13.5 

9 

Highway  49/89 
Sierraville 

6/17/91 

Chip-Seal 

3.5  mi 

Yes 

5 

10 

Highway  49 
Placer  Cnty 

6/8/92 

Saw 

Cutting 

1,C91  -  ' 
1,526  ft 

No 

2 

11 

Highway  50 
El  Dorado  Cnty 

6/9/92  - 
6/10/92 

Saw 
Cutting 

1,380  ft 

No 

5 

12 

Highway  89 
El  Dorado  Cnty 

6/11/92  - 
6/12/92 

Shoulder 
Work 

2,962  - 
3,421  ft 

No 

6 

13 

Highway  10  8 
Tuolome  Cnty 

6/22/92 

Shoulder 
Work 

985  ft 

No 

5.5 

14 

Highway  88 
Alpine  Cnty 

6/23/92 

Striping 

8,828  ft 

Yes 

3.5 

15** 

Highway  50 
El  Dorado  Cnty 

6/10/92 

A-C 

3,760 

Yes 

3.5 

Pilot  study  site;  observations  collected  in  1  direction  only 
'   Validation  site;  delay  data  collected 
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computer  files  were  programmed  to  display  observation  times  with  notations  defining  the  specific 
event  associated  with  each  recorded  time. 

The  two  data  collectors  positioned  themselves  at  opposite  ends  of  the  work  zone  (along 
side  the  flaggers).  Using  the  lap-top  computers,  the  data  collectors  simultaneously  recorded  the 
actual  times  of  relevant  events.  Specifically,  the  real-time  recorded  observations  were: 

1.  The  times  at  which  the  flaggers  initiated  right-of-way  for  the  subject  direction  of  travel. 

2.  The  times  at  which  each  queued  vehicle  actually  "entered"  the  work  zone  (i.e.,  the  times 
the  front  bumper  of  each  vehicle  crossed  an  imaginary  line  at  the  flagger's  station). 

3.  The  times  at  which  the  last  queued  vehicle  "entered"  the  work  zone. 

4.  The  times  of  any  additional  vehicles  entering  the  work  zone  after  the  queue  had 
dissipated. 

5.  The  times  at  which  the  flagger  terminated  right-of-way  for  the  subject  direction  of  travel. 

6.  The  times  at  which  vehicles  traveling  in  the  opposing  direction  exited  the  designated  work 
zone  area. 

An  example  computer  print-out  for  a  single  cycle  is  presented  in  Figure  3-1. 

The  data  collectors  individually  recorded  these  event  times  at  each  end  of  the  work  zone 
over  repeated  cycles.  These  data  provide  virtually  all  relevant  information  for  defining  motorist 
delays   at  the  work  zone. 

•  The  elapsed  time  between  the  exit  of  the  last  vehicle  in  the  opposing  platoon  and  the 
entry  of  the  first  queued  vehicle  in  the  subject  travel  direction  reflects  start-up  lost  time, 
L. 

•  The  elapsed  times  between  successive  queued  vehicle  entries  into  the  work  zone  represent 
discharge  or  saturation  headway  which  is  the  inverse  of  saturation  flow  rate,  S. 

•  The  elapsed  times  between  entry  of  the  last  queued  vehicle  and  the  final  vehicle  entry  in 
the  same  direction  (each  cycle)  represents  the  green  time  extension,  E. 
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13 :55:33 

Flagman  ir.it 

iates  ROW 

13:55:40 

Car  enters 

13:55:44 

Car  enters 

13:55:48 

Car  enters 

13:55:50 

Car  enters 

13:55:53 

Truck  enters 

13:55:58 

Car  enters 

13:56:01 

Car  enters 

13:56:02 

Car  enters 

13:56:05 

Car  enters 

13 :56:07 

Car  enters 

13:56:10 

Last  CAR  in 

Queue 

13:56:15 

Car  enters 

13:56:18 

Flagman  term 

inates  ROW 

13:59:30 

1st  vehicle 

opposing  direction 

14:02:58 

Last  vehicle 

opposing  direction 

14:03:03 

Flagman  init 

iates  ROW 

Figure  3-1 
Example  Data  File 


Given  the  prevailing  "first-in,  first-out"  queueing  discipline  and  the  absence  of  major 
access  points  within  the  work  zone,  the  sequential  work  zone  entry  and  exit  times  define 
work  zone  travel  times  (and  thus  speeds).  Typically,  vehicle  travel  times  were  sampled 
for  first  and  last  vehicles  in  each  platoon. 
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Some  amount  of  error  likely  occurred  when  measuring  travel  times  based  on  vehicle  entry 
and  exit  times.  Such  errors  could  have  occurred  as  a  result  of  vehicles  joining  or  leaving 
platoons  via  the  minor  access  points  within  the  designated  work  zone  area.  However,  even 
under  extreme  scenarios,  measured  errors  would  generally  lead  to  computed  speed  errors 
of  less  than  1  mph.  This  is  well  within  the  range  of  random  speed  fluctuations  naturally 
occurring  within  the  work  zone.  Moreover,  the  measured  sample  travel  times  were 
augmented  with  numerous  "floating  car  runs"  where  travel  times  were  measured  in  a  test 
car  operated  by  the  data  collection  team. 

•  Cycle  lengths  were  measured  as  previously  defined. 

•  As  queues  often  extended  beyond  view  of  the  data  collectors,  vehicle  arrival  times  were 
generally  not  recorded  simultaneously  with  vehicle  discharge  times.  This  means  that 
hourly  traffic  demands  can  only  be  approximated  (rather  than  precisely  determined)  from 
the  resulting  data  sets.  This  approximation  is  not  a  concern  given  that  relevant  parameters 
(i.e.  work  zone  travel  speeds,  saturation  headways  and  lost  times)  will  not  vary  as  a 
function  of  small  changes  in  traffic  demand.  Moreover,  vehicle  arrival  times  were 
separately  sampled  upstream  of  work  zone  entrances  to  confirm  arrival  patterns. 

Data  were  typically  collected  over  90  minute  time  intervals. 


3.2.  Data  Analysis 

For  illustration  purposes,  table  3-2  presents  directional  hourly  demand  rates  over  each  90- 
minute  observation  interval  measured  at  site  numbers  3,5  and  7  -  each  of  which  represents  work 
zones  operating  under  A-C  overlay  activities. 
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Table  3-2  Observed  Traffic  Demand  Rates 


Site  Number 

Observation 
Time 

Demand  (vph) 

Northbound 

Southbound 

Total 

3 

8:00-9:30  am 
10:00-11:30  am 

201 

227 

165 

243 

366 

470 

5 

8:00-9:00  am 

10:00-11:30  am 

12:00-1:30  am 

217 
272 
227 

167 

225 
258 

384 
497 
485 

7 

1:30-3:00  pm 

215 

302 

517 

Table  3-3  presents  statistics  for  observed  cycle  length  for  each  90-minute  interval  (in 
both  directions  of  travel).  The  observed  operations  exhibit  considerable  variation  from  cycle  to 
cycle.  Referring  to  the  values  of  each  coefficient  of  variation  (COV),  cycle  lengths  observed 
over  each  90-minute  interval  vary  from  12  to  41  percent.  These  variations  are  attributable  to 
fluctuations  in  arrival  and  departure  rates  as  well  as  occasional  interruptions  to  traffic  flow  caused 
by  the  work  activity. 

The  variability  in  cycle  length  is  exhibited  not  only  within  sites,  but  between  sites  as  well. 
Although  conditions  prevailing  at  site  7  do  not  differ  dramatically  from  the  other  two  locations, 
the  site's  average  cycle  lengths  exceed  the  other  two  locations  by  a  factor  of  almost  two. 
Evaluation  of  the  data  collected  at  site  7  indicated  that  the  location  operates  with  significantly 
higher  start-up  lost  times  and  a  higher  frequency  of  random  interruptions  (reflected  in  the 
periodic  occurrence  of  large  individual  saturation  headways). 
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This  observed  variability  in  work  zone  operation  (i.e.  cycle  length)  obviates  the  need  for 
a  probabilistic  approach  to  delay  prediction  where  under-saturated  conditions  prevail. 


Table  3-3  Observed  Cycle  Length  Statistics 


Site  No. 

Observation 
Time 

Travel 
Direction 

Cycle  Length 

Mean 
(mins) 

Variance 
(mins2) 

cov 

3 

8:00-9:30  a.m. 

Northbound 
Southbound 

10.6 
10.6 

3.3 
3.1 

0.17 
0.17 

3 

10:00-11:30  a.m. 

Northbound 
Southbound 

11.7 
12.6 

22.6 
15.4 

0.41 
0.31 

5 

8:00-9:30  a.m. 

Northbound 
Southbound 

10.8 
10.8 

1.9 

3.2 

0.13 
0.17 

5 

10:00-11:30  a.m. 

Northbound 
Southbound 

12.3 
12.1 

3.6 

2.6 

0.15 
0.13 

5 

12:00-1:30  p.m. 

Northbound 
Southbound 

10.5 
10.0 

1.5 

1.7 

0.12 
0.13 

7 

1:30-3:00  p.m. 

Northbound 
Southbound 

20.4 
19.7 

40.6 
40.1 

0.31 
0.32 
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4.    SIMULATION  MODEL 


The  accuracy  of  the  proposed  model  was  evaluated  using  a  microscopic,  stochastic 
computer  simulation  model  developed  as  part  of  this  research  project.  This  report  section 
describes  the  simulation  model  and  simulation-based  validation  efforts.  The  final  portion  of  this 
report  section  documents  the  validation  tasks  exploiting  the  available  empirical  data  base. 


4.1  Microscopic  Simulation  Model 

A  schematic  flow  chart  of  the  FORTRAN-based  simulation  model  is  presented  in  Figure 
4-1.  The  microscopic  simulation  model  is  designed  to  estimate  individual  vehicles'  queueing 
delay  at  the  work  zone  entrance  and  does  not  generate  the  delay  created  by  traveling  through  the 
work  zone  at  less  than  desired  speed.  Queueing  delay  (alone)  is  sufficient  for  validation  purposes 
as  work  zone  travel  delay  is  clearly  the  difference  between  actual  and  desired  travel  times. 

Directional  vehicle  demand  rates,  the  physical  length  of  the  lane  closure  and  the  time 
period  to  be  analyzed  are  specified  inputs.  Directional  vehicle  arrival  times,  queue  discharge 
headways  and  work  zone  travel  speeds  of  last  vehicles  in  platoon  are  microscopically  selected 
by  the  model  based  upon  empirically  identified  distributions.  Following  this,  the  simulation  model 
replicates  the  operational  characteristics  occurring  under  one  way  traffic  control. 

•  Following  the  initiation  of  right-of-way  in  a  given  direction,  queued  vehicles  discharge 
into  the  work  zone  at  so-called  "saturation  flow  rate." 

•  Once  the  directional  queue  has  completely  dissipated,  right-of-way  may  be  extended  to 
serve  additional  vehicles  arriving  "shortly  after"  the  queue  dissipation  time.  Specifically, 
right-of-way  is  continued  until  the  elapsed  time  between  consecutively  arriving  vehicles 
exceeds  12.2  seconds.  The  12.2-second  limit  was  empirically  derived  using  discrete  choice 
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INPUT 

Arrival  Headway 

Discharge  Headway 

Last  Vehicle's  Speed 


I  :  Individual  Vehicles 
K  :  Counting  Number  of 
Vehicles  in  each  Cycle 


Set  Clock  =  0.0 


Time  ROW  Ends 


1=1 


K=l 


Calculate 
Departure  Time 


No 


Calculate 
Queueing  Delay 


1=1+1,  K=K+1 


Departure  Time 
=  Arrival  Time 


ROW  End  Time  =  Departure  time  of  Last 
Vehicle  of  The  Cycle 
Compute  Average  Delay  per  Cycle 


Figure  4-1 
Microscopic  Simulation  Model 
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analysis  [Ben-Akiva  and  Lerman,  1985] 

All  discharged  vehicles  traverse  the  work  zone  segment. 

Once  the  "last  vehicle"  in  the  discharging  direction  has  completely  traversed  the  work 
zone,  some  minor  "lost  time,"  on  the  order  of  a  few  seconds,  occurs  prior  to  right-of-way 
allocation  to  the  opposing  directional  queue.  (Given  it's  small  contribution  to  operating 
time,  this  lost  time  was  deterministically  modeled  in  the  simulation).  The  process  is  then 
repeated  for  the  opposing  directional  movement.  Several  hours  of  simulation  time  were 
used  to  initialize  (i.e.,  "warm-up")  each  experiment  to  insure  the  system  (i.e.,  work  zone) 
had  reached  steady  state  conditions  [Son,  Cassidy  &  Madanat,  1994]. 


4.1.1.  Threshold  Headways 

During  one-way  traffic  control,  the  extension  of  right-of-way  beyond  queue  dissipation 
reflects  the  subjective  judgement  of  flaggers.  This  decision  process  can  be  classified  as  a  binary 
choice  to  extend  or  terminate  green  time  based  upon  the  headway  of  the  next  arriving  vehicle. 
Thus,  the  maximum  headway  that  motivates  the  flagger  to  extend  green  time  (i.e.,  the  threshold 
headway)  was  identified  and  introduced  in  the  simulation  model. 

A  probit  model  [Ben-Akiva  and  Lerman,  1985]  was  likewise  estimated  from  empirical 
observations  to  identify  the  probability  of  right-of-way  termination  (by  flaggers)  following  queue 
dissipation  as  a  function  of  the  next  arriving  vehicle's  headway  length.  A  mean  value  of  H  (a 
"threshold"  headway)  was  derived  from  the  probit  function  [Mahmassani  and  Sheffi,  1981]. 

The  model  recognizes  variations  in  threshold  values  across  flaggers.  All  headways  which 
are  larger  than  the  flagger' s  threshold  headway  are  invariably  rejected  (right  of  way  terminated). 
Conversely,  flaggers  accept  all  headways  smaller  than  or  equal  to  the  threshold  headway.  It  is 
further  assumed  that  these  threshold  values  follow  some  probability  distribution  across  the 
population  of  flaggers.  The  mean  of  this  distribution  is  the  relevant  parameter  in  this  model. 
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If  we  define  Pr(a) ,  as  the  probability  that  a  headway  t  after  queue  dissipation  is  greater 
than  or  equal  to  the  flagger's  threshold  value  Hlh, 


Pr{a)  =  P{t  >  HJ 

=  Pit  >  Hth+e) 
=  Pie  <  t-HJ 


(4.1) 


Where  Hlh  is  the  mean  value  of  thresholds  and  e  is  the  random  deviation  or  disturbance  of  each 

flagger's  threshold  value  from  the  mean.  One  logical  assumption  is  to  view  the  disturbance  as 
the  sum  of  a  large  number  of  unobserved  and  independent  components.  By  the  central  limit 
theorem,  the  distribution  of  the  disturbances  would  tend  to  be  normal. 


Piia)  =  prie  <  t-H,) 


=  O 


f     -  \ 
t-H., 


v    a    J 


(4.2) 


where  a  is  the  standard  deviation  of  the  disturbance  and  <£(•)  denotes  the  standardized  cumulative 
normal  distribution. 

Equation  (4.2)  can  be  rewritten  as 


where: 


«     ■    P#*. 


Prid)  =  Qm-HJ] 
=  <D[pr-a] 


(4.3) 
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Values  of  6   and   3    (tne  estimates  of  a   and  (3)  are  obtained  through  Maximum 
Likelihood  Estimation  [Benjamin  and  Cornell,  1970].  The  value  of  Hth  is  estimated  using  the 


relation 


m  -& 
,h'J 


(4.4) 


Right-of-way  is  invariably  extended  whenever  an  arriving  vehicle's  headway  is  less  than  Hth  and 
right-of-way  is  terminated  when  a  headway  exceeds  Hth.  Figure  4-2  presents  the  Maximum 
Likelihood  Estimation  results.  Using  Equation  (4.4),    H h  is  estimated  to  be  12.2  seconds. 


Independent 

Estimated 

Standard 

t  -  Statistics 

Variable 

Coefficient 

Error 

one 

2.32083 

0.25992 

8.92911 

headway 

-0.19006 

2.62013e-002 

-7.25379 

Auxiliary  Statistics 

at  Convergence 

Initial 

log  likelihood 

-79.622 

-155.96 

number  of  observation 

225 

percent  correctly  predicted 

82.667 

Figure  4-2 
Maximum  Likelihood  Estimation 
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4.2.  Evaluation  of  Model  Accuracy 

As  a  means  of  evaluating  the  closed-form  expressions  presented  in  section  2  of  this  report, 
delay  estimates  yielded  from  these  models  were  compared  with  simulated  values  for  a  range  of 
specified  operating  conditions. 


4.2.1.  Experimental  Design 

Eighteen  work  zone  operating  scenarios  were  selected  to  represent  commonly  observed 
conditions  in  the  field.  The  18  scenarios,  presented  in  Table  4-1,  reflect  a  range  of  directional 
traffic  demand  rates  and  work  zone  physical  lengths.  Each  of  the  18  scenarios  were  simulated 
for  1,000  operating  cycles  (following  a  suitable  "warm-up"  timej  to  generate  average  values  of 
directional  delay  at  the  work  zone  entrances.  As  the  microscopic,  stochastic  simulation  model 
replicates  the  actual  phenomena  of  work  zone  operation  and  all  relevant  random  variables 
conform  to  empirically  observed  distributions,  outcomes  from  the  microscopic  simulations  can 
be  viewed  as  "exact"  (i.e.,  real-world)  solutions. 

Delay  values  generated  through  simulation  were  compared  with  estimates  from  the  closed- 
form,  stochastic  models  in  section  2.  As  green  time  extensions,  Et  represent  relatively  small 
portions  of  the  cycle  length,  delay  estimates  were  derived  from  the  closed-form  expression  by 
1)  incorporating  estimated  values  of  E.  as  described  in  section  2.4.2,  and  2)  assuming  E.  =0. 


4.2.2.  Validation  Findings 

Table  4-2  presents  estimates  of  average  directional  delay  at  the  work  zone  entrance  for 
each  of  the  18  scenarios.   Referring  to  Table  4-2,  differences  between  delay  estimates  generated 
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Table  4-1 


Scenarios  Evaluated  Through  Simulation 


Scenario 
No. 

Work  Zone 

Section  Length 

(miles) 

Flow  Rate 

Direction  1 

(vph) 

Flow  Rate 

Direction  2 

(vph) 

Total  Flow 
Rate 
(vph) 

1 

100 

100 

200 

2 

100 

200 

300 

3 

100 

300 

400 

4 

300 

100 

400 

5 

200 

300 

500 

6 

250 

250 

500 

7 

300 

200 

500 

8 

300 

300 

600 

9 

1.5 

100 

100 

200 

10 

1.5 

100 

200 

300 

11 

1.5 

100 

300 

400 

12 

1.5 

200 

300 

500 

13 

1.5 

300 

300 

600 

14 

0.75 

100 

100 

200 

15 

0.75 

100 

200 

300 

16 

0.75 

100 

300 

400 

17 

0.75 

200 

300 

500 

18 

0.75 

300 

300 

600 
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Table  4-2 


Estimated  Average  Delays 


Scenario 

No. 

Simulation 
(sec) 

(1) 
Stochastic  Model 

Incorporating  E; 

(sec) 

(2) 

Stochastic  Model 

E,=0 

(sec) 

Difference  between 

(1)  and  (2) 

(%) 

1 

279 

317 

319 

0.6 

2 

319 

342 

342 

0.0 

3 

360 

378 

375 

0.8 

4 

300 

342 

346 

1.1 

5 

383 

400 

397 

0.8 

6 

377 

387 

386 

0.3 

7 

348 

376 

376 

0.0 

8 

415 

429 

426 

0.7 

9 

412 

459 

462 

0.6 

10 

456 

491 

492 

0.2 

11 

528 

538 

536 

0.3 

12 

561 

568 

566 

0.4 

13 

600 

608 

605 

0.5 

14 

224 

247 

248 

0.4 

15 

257 

268 

267 

0.4 

16 

283 

298 

295 

1.0 

17 

311 

316 

313 

0.9 

18 

372 

340 

336 

1.0 
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by  simulation  and  by  the  stochastic  models  are  generally  less  than  10  percent. 

Invariably,  the  differences  in  delay  as  estimated  with  and  without  consideration  of    E.  are 

negligible.  To  some  extent,  the  insignificance  of  E.  might  be  attributable  to  pilot  car  operation. 

Where  a  pilot  car  is  utilized,  platoon  speeds  through  the  work  zone  are  generally  much  lower 
than  in  the  absence  of  a  pilot  car  [California,  1978;  Cassidy,  Son  and  Rosowsky,  1993]  as 
platooned  motorists  are  constrained  to  travel  at  speeds  selected  by  the  pilot  car  operator. 
Vehicles  entering  the  work  zone  after  queue  dissipation  typically  travel  at  "desired"  speeds  until 
joining  the  tail  of  the  slower-moving  platoon.  As  these  "late  arrivals"  subsequently  select 
relatively  short  headways,  their  presence  creates  only  minor  additional  time  to  clear  the  work 
zone  before  allocating  right-of-way  to  the  opposing  direction. 

In  the  absence  of  pilot  car  operation,  however,  vehicles  allowed  entry  into  the  work  zone 
following  queue  dissipation  often  do  not  "catch  up"  with  the  through-moving  platoon.  Thus,  the 
initiation  of  right-of-way  for  the  opposing  direction  is  postponed  (until  the  final  discharging 
vehicle  completely  traverses  the  work  zone)  which  may  substantially  increase  delay. 


4.2.3.  Empirical  Validation  Data 

Although  the  microscopic  simulation  model  closely  replicates  "real-world"  conditions, 
empirical  data  perhaps  facilitate  the  most  convincing  validations.  The  limitation  in  applying 
empirical  data  for  the  validation  of  predicted  delay  distributions  stems  from  an  inability  to  collect 
a  sufficient  sample  size.  A  sufficient  data  collection  is  not  possible  given  that  prevailing  work 
zone  operating  conditions  (specifically  vehicle  arrival  rates)  typically  vary  from  hour  to  hour. 
Thus,  conditions  can  not  be  fixed  for  a  time  period  required  to  observe  100  or  more  cycles.  In 
short,  a  reasonable  distribution  of  average  cyclic  delays  can  not  be  observed  in  the  field. 

One  can,  however,  field  measure  directional  delays  for  finite  periods  of  time  and  compare 
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empirically  observed  average  values  with  average  values  generated  form  the  queueing  models. 
Such  comparisons  would  verify  the  predictive  capabilities  of  the  models  which  lie  at  the 
foundation  of  the  methodologies  for  both  under-saturated  and  over-saturated  conditions. 
The  site  used  for  collecting  validation  data  was  Highway  50  in  El  Dorado  county,  California. 
Approximately  90  minutes  of  vehicle  delay  observations  were  collected  in  each  direction  of 
travel.  Individual  vehicle  delays  were  measured  as  the  difference  between  actual  and  desired 
travel  times  between  the  work  zone  entrance  and  a  fixed  location  upstream  of  queueing  activity. 
Computed  free-flow  travel  times  were  subtracted  from  the  actual  measured  travel  times 
approaching  the  work  zone  entrance.  Average  values  were  reflected  as  the  summed  total  delay 
values  divided  by  the  number  of  arriving  vehicles. 

The  field-measured  values  of  average  delay  in  queue  for  each  direction  of  travel  are 
compared  with  deterministically  predicted  values  in  Table  4-3.  Discrepancies  are  under  10 
percent. 


Table  4-3 
Field-Measured  and  Predicted  Average  Delay  with  Queueing  Model  (for  subject  direction) 


Length 
(ft) 

Traffic  Volume  (vph) 

Field- 
Measured 
(mins/veh) 

Predicted 
(mins/veh) 

Error 
(%) 

subject  dir. 

opposite  dir. 

3,760 

309 

207 

3.7 

4.1 

9.7 

3,760 

396 

310 

5.0 

5.5 

9.0 

3,760 

328 

390 

5.9 

6.3 

6.3 
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5.  SUMMARY  AND  CONCLUSIONS 


This  report  has  described  tasks  associated  with  the  development  of  a  methodology  for 
assessing  traffic  flow  impacts  (i.e.,  delay)  resulting  from  the  implementation  of  one-way  traffic 
control.  Methodologies  have  been  formulated  for  both  under-saturated  and  over-saturated 
operating  conditions. 

Stochastic  queueing  models  are  used  to  estimate  average  delay  prevailing  at  two-lane 
highway  work  zone  for  under-saturated  conditions.  The  models  are  based  upon  previously  derived 
expressions  for  delay  estimation  at  the  V-A  signalized  intersection  of  one-way  streets  [Newell, 
1969]  adapted  to  accommodate  operating  characteristics  unique  to  one-way  traffic  control.  The 
models  capture  stochastic  affects  on  expected  delay  as  a  function  of  directional  traffic  demands, 
work  zone  physical  length  and  observable  traffic  stream  measures  prevailing  at  under-saturated 
work  zone  operations. 

For  over-saturated  conditions,  increased  delays  created  by  random  effects  (e.g.,  variations 
in  traffic  stream  characteristics)  are  very  small  relative  to  the  deterministic  delay  components. 
As  such,  the  proposed  methodology  for  assessing  over-saturated  operation  exploits  deterministic 
queueing  analysis  techniques. 

The  models  presented  in  this  report  estimate  work  zone  queueing  delay.  The  additional 
delay  which  is  incurred  by  motorists  as  they  travel  through  the  work  zone  at  less  than  desired 
speed  is  merely  the  difference  between  the  expected  work  zone  travel  times  and  a  presumed 
travel  time  in  the  absence  of  a  work  zone. 

To  evaluate  the  accuracy  of  the  proposed  stochastic  models,  a  computer  simulation  model 
was  developed.  Based  upon  the  outcomes  generated  by  the  simulation  model,  the  proposed  delay 
expressions  appear  to  adequately  reflect  the  impacts  of  the  work  zone  operation  at  a  two-lane 
highway.  Validation  tasks  using  empirical  data  confirm  this  conclusion. 

Finally,  to  simplify  and  expedite  the  application  of  the  proposed  methodologies,  ah 
formulated  procedures  have  been  fully  computerized.  The  resulting  software  package  facilitates 
application  of  the  methodologies  with  a  minimum  of  time  and  effort.  A  software  User's  Manual 
is  contained  in  Appendix  B  of  this  report. 
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APPENDIX  A:  PARAMETER  VALUES  FOR  EACH  WORK  ACTIVITY  TYPE 

Table  A-l  presents  field-measured  values  for  relevant  operational  parameters  observed 
under  all  four  types  of  work  activity.  The  upper  portion  of  this  table  presents  average  parameter 
values  and  the  lower  portion  tabulates  observed  standard  deviations  required  for  estimating  the 
average  delay  per  cycle. 


Table  A-l 
Model  Parameter  Values 


Parameters 

WORK  ACTIVITY  TYPE 

A-C 

Chip-Seal 

Construction-a 

Construction-b 

Expected  value 

Saturation  Flow 
(vph) 

1084 

1018 

1090 

1048 

Speed 
(ft/sec) 

34 

29 

45 

30 

Lost  Time 
(sees) 

30 

23 

35 

15 

Variance 

Speed 
(ft/sec)2 

10.52 

7.72 

11.02 

3.22 

Lost  Time 

(sees)2 

0 

0 

104.92 

0 

Threshold 
Headway(secs) 

12.21 

—  . 

a-  without  pilot  car 
b-  with  pilot  car 
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APPENDIX  B:  USER  MANUAL 
work  -  Software  for  Analyzing  One-way  Traffic  Control  Impacts 


B.l  Introduction 

The  computer  program  "work"  has  been  developed  to  assess  the  impacts  imparted  to 
vehicles  as  a  result  of  one-way  traffic  control  implementation.  The  program,  which  can  evaluate 
both  undersaturated  and  oversaturated  operating  conditions,  exploits  the  computational  procedures 
previously  described  in  the  final  report.  The  reader  should  refer  to  the  report  for  procedural 
details  concerning  the  delay  prediction  models  as  well  as  any  technical  issues  noted  in  this 
Appendix  (i.e.,  User  Manual).  As  the  name  implies,  the  subject  appendix  addresses  use  of  the 
software  packages  itself. 


B.2  Installing  the  work  Software 


The  Software  can  be  run  on  IBM-compatible  personal  computers.  A  math  co-processor 
is  required  where  work  is  to  be  operated  on  286-machines  (i.e.,  80286). 

Although  the  work  software  can  be  assessed  and  executed  from  the  diskette,  computational 
speed  is  enhanced  by  accessing  the  program  from  the  "hard  disk".  Copying  the  contents  of  the 
work  diskette  to  the  hard  disk  can  be  performed  using  standard  DOS  commands.  The  user  may 
wish  to  first  create  a  unique  directory  using  the  commands 

mtbdirname 
where  dirname  represents  the  user-specified  name  of  the  directory.  Once  changing  to  this 
directory  using  the  command 

cd\dirname 
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the  contents  of  the  work  diskette  can  be  transferred  to  the  specified  directory  with  the  command 

copy  rf:work.exe 
where  d  represents  the  drive  where  the  work  diskette  has  been  inserted. 


B.3  Executing  work  -  Program  Input 

The  program  is  accessed  by  typing  work.  Once  accessing  the  program,  the  user  is  merely 
required  to  respond  to  all  prompts  by  typing  in  the  appropriate  information  and  then  hitting 
<return/enter>: 

a)  name  of  highway  (up  to  10  characters) 

This  information  is  optional,  although  such  data  can  aid  in  "book-keeping"  matters. 

b)  if  the  highway  is  east-west,  type  1 
if  the  highway  is  north-south  type  2 

This  information  is  used  to  "orientate"  the  analysis  to  the  directional  alignment  of  the 
subject  highway.  Typing  a  numeric  character  other  than  1  or  2  will  not  effect  procedural 
computations.  However,  entering  a  number  other  than  1  or  2  will  mean  that  subsequent  prompts 
(and  print-out)  will  not  display  travel  directions. 

c)  hours  of  analysis 

This  prompt  represents  the  total  number  of  consecutive  hours  over  which  analysis  is  to 
occur.  Selecting  a  number  for  hours  of  analysis  is  not  particularly  important  where  operation 
reflects  undersaturated  conditions  given  that,  under  these  conditions,  operation  during  a  previous 
hour  will  not  influence  the  following  hour.  However,  where  oversaturated  conditions  prevail,  the 
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occurrence  of  residual  queueing  impacts  operation  from  one  hour  to  the  next.  The  reader  may 
refer  to  "Example  2"  in  this  appendix  for  a  brief  discussion  of  this  issue. 

Up  to  29  consecutive  hours  can  be  specified  for  evaluation  using  work.  Specifying  more 
than  29  hours  will  eventually  result  in  a  "run  time  error". 

d)  maximum  green  time  for  northbound  or  eastbound  direction  (mins) 

This  information  reflects  the  assumed  maximum  time  (in  minutes)  over  which  flaggers 
are  willing  to  allocate  right-of-way  before  "cutting-off '  a  discharging  queue  of  vehicles.  This 
specified  value,  which  should  be  based  upon  judgement  or  existing  operating  policies,  is 
important  only  in  the  context  of  oversaturated  operating  conditions.  Where  a  prevailing 
directional  green  time  calculated  by  the  software  exceeds  the  maximum  specified  value, 
oversaturated  conditions  will  occur.  The  directional  green  times  assumed  to  prevail  during 
oversaturated  conditions  are  the  specified  maximum  values. 

Specified  maximum  green  time  values  in  work  must  be  greater  than  0  (minutes)  and  less 
than  60  (minutes). 

e)  maximum  green  time  for  southbound  or  westbound  direction  (mins) 
same  as  above 

f)  work  zone  length  (miles) 

The  physical  length  of  the  work  zone  (in  miles)  is  to  be  entered  here.  Any  value  of  length 
can  be  specified  where  work  activity  is  A-C,  chip-seal  or  general  construction  (i.e.,  "other")  with 
the  use  of  a  pilot  car.  The  somewhat  unique  operating  characteristics  of  general  construction 
activity  without  the  use  of  a  pilot  car  dictates  that  specified  work  zone  lengths  should  be  greater 
than  0.25  (miles)  for  avoiding  erroneous  predictions. 

g)  vehicle  demand  rate  -  northbound  or  eastbound  (veh/hr) 
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hour  1 


The  user  is  required  to  enter  directional  vehicle  demand  rate  for  each  hour  of  analysis. 
The  user  is  sequentially  prompted  for  hourly  demand  rate  up  to  the  hour  number  previously 
specified  in  c).  The  prompts  sequentially  request  all  hourly  demand  rate  entries  for  the 
northbound  or  eastbound  directions  before  requesting  the  equivalent  information  for  the  opposing 
travel  direction. 

work  can  accommodate  any  specified  demand  rate. 

h)         speed  limit  or  estimated  prevailing  speed  in  the  absence  of  a  work  zone  (mile/hr) 

An  estimate  of  the  subject  highway's  "normal"  prevailing  average  speed  in  the  absence 
of  one-way  traffic  control  is  required  to  determine  travel  time  delay  through  the  work  zone  itself. 

work  will  accommodate  any  specified  value  of  "normal"  speed.  However,  specified  values 
less  than  35  (mph)  may  result  in  negative  travel  time  delays  as  observed  work  zone  vehicle 
speeds  for  some  job  types  approached  35  mph. 

i)         job  type 

asphalt  concrete  =  1 

chip  seal  =  2 

other  if  pilot  car  is  not  used  =  3 

other  if  pilot  car  is  used  =  4 

The  user  is  required  to  specify  the  type  of  work  activity  prevailing  at  the  work  zone  by 
selecting  options  (i.e.,  keys)  1,2,3  or  4.  This  is  important  input  information  as  parameter  values 
such  as  vehicle  queue  discharge  rates  and  work  zone  travel  speeds  will  vary  as  a  function  of  job 
type. 

The  term  "other"  refers  to  general  construction  activity  such  as  shoulder  work,  cross-drain 
installation  or  restriping  activity.  Where  a  job  type  to  be  analyzed  does  not  fall  into  one  of  the 
four  available  categories,  the  user  should  select  the  job  type  option  most  closely  conforming  to 
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the  work  activity  in  question. 

j)  Comment  (up  to  30  characters) 

This  prompt  provides  the  user  with  the  option  to  input  relevant  characteristics  of  the 
subject  analysis  as  a  means  to  assist  in  "book-keeping".  The  user  can  enter  up  to  30  characters 
(including  spaces  between  characters.) 

k)         File  name  missing  or  blank  -  please  enter  file  name 

This  prompts  is  merely  requesting  the  user  to  name  the  output  file  which  will  represent 
the  analysis  of  the  now-specified  input  conditions.  Filenames  can  consist  of  up  to  8  alpha  and/or 
numeric  characters. 

1)  Do  you  want  to  change  options  ?  (Yes  =  y,  No  =  n) 

This  prompt  represents  the  work  program's  only  editing  function.  Specifically,  the  prompt 
is  intended  to  facilitate  sensitivity  analyses.  By  responding  "y"  to  this  prompt,  the  user  can 
change  one  or  more  operating  parameters  (and/or  comments)  and  thereby  direct  the  work  program 
to  automatically  re-analyze  operation  under  modified  scenarios. 

Parameters  which  can  be  directly  modified  through  this  prompt  are 

1.  work  zone  length 

2.  maximum  directional  green  times 

3.  "normal"  speed  in  the  absence  of  a  work  zone 

4.  job  type 

The  benefits  offered  by  this  "editing  function"  are  perhaps  obvious.  Assume,  for  example, 
that  a  user  is  interested  in  identifying  a  work  zone  physical  length  which  will  result  in 
"acceptable"  delays.  Assume  also  that  the  user  wishes  to  evaluate  operation  over  numerous 
consecutive  hours.    Rather  than   re-entering    all   required   input  information  (including  hourly 
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demands)  and  re-executing  the  work  program  under  varying  section  lengths,  the  user  can  a  priori 
specify  an  array  of  possible  work  zone  lengths  for  the  previously  specified  input  conditions  and 
automatically  execute  an  analysis  for  each  specified  length. 

The  user  can  specify  desired  input  modifications  by  appropriately  responding  to  each 
prompt.  Once  the  user  has  "cycled  through"  one  sequence  of  prompts,  the  software  will  again 
display  the  prompt  Do  you  want  change  options  ?  (  Yes  =  y,  No  =  n  ).  In  this  way,  the  user 
can  perform  virtually  unlimited  iterations  under  modified  input  conditions.  Each  sequence  of 
prompts  will  automatically  result  in  a  separate  analysis  reflecting  specified  changes. 

Figure  B-l  presents  the  computer  screen  reflecting  the  input  information  required  for 
Example  Problem  1. 


B.4  Program  Output 

Promptly  after  responding  to  the  editing  feature  (and  hitting  <return>),  the  work  program 
will  execute  all  computations.  The  resulting  message  display  is 

Stop  -  Program  terminated 

The  output  file  (i.e.,  the  resulting  analyses)  can  be  viewed  directly  from  the  computer  monitor 
and/or  permanent  output  can  be  printed  using  a  standard  printer  device. 

To  view  output  directly  on  the  computer,  the  user  can  type  the  DOS  command 

type  filename 

where  filename  represents  the  name  assigned  to  the  output  file  described  in  k).  This  command 
will  cause  the  entire  output  to  scroll  past  the  monitor,  stopping  at  the  end  of  the  file. 
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name  of  highway   (up  to  10  characters) 

Hwy  4  9  Down 

if  the  highway  is  east-west,  type  1 

if  the  highway  is  north-south,  type  2 

2 

hours  of  analysis 

2 

maximum  green  time  for  northbound  direction  (mins) 

7 

maximum  green  time  for  southbound  direction  (mins) 

7 

work  zone  length  (miles) 

3.25 

vehicle  demand  rate-northbound  (veh/hr) 

hour        1 

170 

hour        2 

160 

vehicle  demand  rate- southbound  (veh/hr) 

hour        1 

140 

hour        2 

145 

speed  limit  or  estimated  prevailing  speed 

in  the  absence  of  a  work  zone  (mile/hr) 

45 

job  type 

Asphalt  Concrete  Overlay  =  1 

Chip  Seal  =  2 

other  if  pilot  car  is  not  used      =  3 

other  if  pilot  car  is  used  =  4 

2 

Comment  (up  to  3  0  characters) 

Example  Problem  1 

File  name  missing  or  blank  -  please  enter  file  name 

UNIT  2?  Testl 

Do  you  want  to  change  options  ?  (Yes  =  y  ,  No  =  n  ) 

n 

Stop  -  Program  terminated 


Figure  B-l 
Example  Input  Information 
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To  review  the  output  file  on  the  monitor  page  by  page  requires  the  DOS  command 

type  filename\more 
Each  new  output  page  is  displayed  on  the  monitor  by  hitting  <return>. 
Printing  a  "hard  copy"  of  the  output  requires  the  DOS  command 

print  filename 

Where  the  printer  has  not  been  previously  used  to  print  work  output,  the  analyst  is  required  to 
initialize  the  printer  by  hitting  <return>  in  response  to  the  prompt 

name  of  list  device  [prn] 

The  prompt  filename  is  currently  being  printed  is  displayed  whenever  a  print-out  is 
requested  by  the  user. 


B.5  Sample  Output  /  Example  Problems 

The  format  and  the  information  conveyed  on  work  output  depend  upon  whether  prevailing 
operation  is  undersaturated  or  oversaturated.  These  format  differences  are  consistant  with  the 
differences  in  actual  computations.  The  differences  are  illustrated  in  three  Example  Problems. 

Example  Problem  1  :  undersaturated  conditions 
Required  Input  Information 
Location  :  Highway  49,  Downeville 
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Orientation  :  North-South 

hours  of  analysis  :  2  (e.g.  8:00  am  through  10:00  am) 

maximum  green  time  for  both  directions  of  travel  :  7  minutes 

work  zone  length  :  3.25  miles 

demand  rates  (vph)  : 

Northbound     Southbound 
Hour  1  170  140 

Hour  2  160  145 

Speed  limit  / prevailing  "normal"  speed  :  45  mph 
job  type  :  chip  seal 

The  input  information  (displayed  on  the  computer  monitor)  for  Example  Problem  1  has 
been  previously  presented  in  Figure  B-l.  Figure  B-2  presents  the  output  generated  form  the  work 
program.  Referring  to  Figure  B-2,  relevant  input  information  is  reproduced  on  the  output  file. 
Both  hours  of  analysis  are  undersaturated.  For  each  of  the  two  hours,  the  average  delay  per 
vehicle  (in  minutes)  and  total  delay  (in  vehicle-hours)  is  displayed  for  each  travel  direction. 

If  the  user  so  desires,  total  directional  delay  can  be  calculated  as  the  sum  of  each 
directional  hourly  value.  Total  average  delay  in  both  directions  is  a  weighted  average  of  hourly 
directional  delay  and  hourly  directional  demand  rates.  The  output  also  displays  initial  and  residual 
queue  lengths  (measured  in  vehicles)  which  are  zero  for  undersaturated  conditions. 

ANALYSIS  RESULTS 


Highway 

Hours  of  Analysis 

Job  Type 

Work  Zone  Length 

Maximum  Green  Time 
northbound 
southbound 


Hwy 4  9  Down 
2     hour(s) 
Chip  Seal 

3.25    mile(s) 


7.00    min 
7.00    min 


Figure  B-2 
Output  -  Example  Problem  1 
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Normal  Speed 
Comments    : 


45.00    mph 
Example  Problem  1 


Vehicle  Demand  Rate  (veh/hr) 


northbound 

southbound 

hour  1 

170. 

140. 

hour  2 

160. 

145. 

hour  1   under- saturated 


northbound 
southbound 


Initial  Queue 

0. 
0. 


Residual  Queue 

0. 
0. 


Delay   -   Average  (min/veh) 


northbound 
southbound 


19.1 
19.5 


Total  (veh-hr) 

54.05 
45.38 


hour  2 


under- saturated 


northbound 
southbound 


Initial  Queue 

0. 
0. 


Residual  Queue 

0. 
0. 


Delay   - 

northbound 
southbound 


Average  (min/veh) 

19.1 
19.3 


Total  (veh-hr) 

51.04 
46.57 


hour  1 
hour  2 


Total  Delay  (veh-hrs) 


northbound 
54.05 
51.04 


southbound 
45.38 
46.57 


Figure  B-2  (con't) 
Output  -  Example  Problem  1 
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Example  Problem  2  :  oversaturated  conditions 

Required  Input  Information 

Location  :  Highway  28,  Tahoe  City 

Orientation  :  East-West 

hours  of  analysis  :  3  (e.g.  2:00  pm  through  5:00  pm) 

maximum  green  time  for  both  directions  of  travel  :  5  minutes 

work  zone  length  :  1.25  miles 

demand  rates  (vph)  : 

Eastbound        Westbound 
Hour  1  440  355 

Hour  2  460  375 

Hour  3  490  360 

speed  limit  / prevailing  "normal"  speed  :  50  mph 
job  type  :  A-C 

ANALYSIS  RESULTS 


Highway 

Hwy 

28 

Hours  of  Analysis 

3    : 

hour ( s ) 

Job  Type 

Asphalt  Concrete 

Overlay 

Work  Zone  Length 

1.25    ] 

nile (s) 

Maximum  Green  Time 

eastbound 

5.00 

min 

westbound 

5.00 

min 

Normal  Speed 

50.00 

mph 

Comments    : 

Example  Problem  2 

Vehicle  Demand  Rate  (veh/hr) 


eastbound 

westbound 

hour  1 

440. 

355. 

hour  2 

460. 

375. 
Figure  B-3 

Output  -  Example  Problem  2 
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hour  3  490.  360. 


hour  1   over-saturated 


Initial  Queue      Residual  Queue 

eastbound  0.  131. 

westbound  0.  46. 


Delay   -   Average  (min/veh)     Total  (veh-hr) 

eastbound  14.6  106.73 

westbound  10.9  64.23 


hour  2   over-saturated 


Initial  Queue      Residual  Queue 

eastbound  131.  282. 

westbound  46.  112. 


Delay   -   Average  (min/veh)     Total  (veh-hr) 

eastbound  25.1  247.51 

westbound  17.1  120.01 


hour  3   over-saturated 


Initial  Queue      Residual  Queue 

eastbound  282.  462. 

westbound  112.  162. 


Delay   -   Average  (min/veh)     Total  (veh-hr) 

eastbound  32.1  413.29 

westbound  22.7  178.29 


Figure  B-3  (con't) 
Output  -  Example  Problem  2 
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Total    Delay    (veh-hrs) 


eastbound 

westbound 

hour  1 

106.73 

64.23 

hour  2 

247.51 

120.01 

hour  3 

413.29 

178.29 

Figure  B-3  (con't) 
Output  -  Example  Problem  2 


As  illustrated  in  Figure  B-3,  each  of  the  three  analysis  hours  in  Example  Problem  2 
operate  at  oversaturated  conditions.  To  further  characterize  operating  conditions  during 
oversaturated  conditions,  the  output  displays  initial  and  residual  queue  lengths  (measured  in 
vehicles)  prevailing  at  the  beginning  and  end  of  each  analysis  hour. 

The  inclusion  of  the  "initial  queue"  at  each  analysis  hour  serves  to  illustrate  a  rather 
important  point.  Computations  assume  that  hours  directly  preceding  selected  analysis  hours  are 
characterized  by  undersaturated  conditions.  Where  this  assumption  is  not  true,  significant 
discrepancies  may  exist  between  actual  and  computed  delay  values.  To  eliminate  these  errors, 
specified  hours  of  analysis  should  include  hours  at  or  before  the  onset  of  oversaturated 
conditions. 


Example  Problem  3  :  under-  and  oversaturated  conditions 

Required  Input  Information 

Location  :  Highway  49,  Placer  County 

Orientation  :  North-South 

hours  of  analysis  :  2  (e.g.  3:00  pm  through  5:00) 
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maximum  green  time  for  both  directions  of  travel  :  5  minutes 
work  zone  length  :  0.75  miles 
demand  rates  (vph)  : 


Eastbound 

Westbound 

Hour  1      260 

210 

Hour  2     495 

330 

speed  limit  /prevailing  "normal"  speed  :  50  mph 
job  type  :  with  pilot  car 

Referring  to  the  output  illustrated  in  Figure  B-4,  the  first  analysis  hour  in  Example 
Problem  3  results  in  undersaturated  operation.  However,  a  surge  in  demand  rate  during  the 
second  hour  produces  oversaturated  conditions  for  that  hour.  Referring  to  Figure  B-4,  a  rather 
large  residual  queue  of  138  vehicles  prevails  in  the  northbound  direction.  Thus,  oversaturated 
conditions  will  be  "carried-over"  into  the  third  hour  and  perhaps  beyond. 

ANALYSIS  RESULTS 


Highway 

Hwy  49 

Hours  of  Analysis 

2 

hour (s) 

Job  Type 

Other   with  Pi 

lot  Car 

Work  Zone  Length 

.75 

mile(s) 

Maximum  Green  Time 

northbound 

5.00 

min 

southbound 

5.00 

min 

Normal  Speed 

50.00 

mph 

Comments    : 

Example  Problem  3 
Figure  B-4 

Output  -  Example  Problem  3 
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Vehicle  Demand  Rate  (veh/hr) 


northbound 

southbound 

hour  1 

260. 

210. 

hour  2 

495. 

330. 

hour  1 


under -saturated 


northbound 
southbound 


Initial  Queue 

0. 
0. 


Residual  Queue 

0. 
0. 


Delay   - 

northbound 
southbound 


Average  (min/veh) 

5.0 
5.2 


Total  (veh-hr) 

21.66 
18.37 


hour  2   over-saturated 


northbound 
southbound 


Initial  Queue 

0. 
0. 


Residual  Queue 

138. 
0. 


Delay   -   Average  (min/veh) 


northbound 
southbound 


12.8 
6.2 


Total  (veh-hr) 

105.47 
34.36 


Total  Delay  (veh-hrs) 


hour  1 
hour  2 


northbound 

21.66 

105.47 


southbound 
18.37 
34.36 


Figure  B-4 
Output  -  Example  Problem  3 
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B.6  Final  Comments  -  User  Notes 

The  work  software  program  was  developed  with  virtually  one  objective  in  mind:  To 
provide  an  accurate  and  easy-to-use  tool  for  evaluating  one-way  traffic  control  scenarios.  To 
insure  that  the  software  would  enjoy  broad-based  implementation,  work  was  designed  to  be 
straightforward  to  use.  By  minimizing  execution  tasks,  the  software  can  be  operated  with  little 
or  no  training  requirements  for  the  user.  As  such,  work  represents  a  very  simple  software  design. 

•  The  program  facilitates  only  limited  editing  features.  (Specifically,  the  capability  to 
perform  sensitivity  analyses  by  responding  yes  to  the  prompt  "  Do  you  want  to  change 
options?")  However,  the  user  is  not  able  to  move  "back  and  forth"  within  the  input  file 
once  data  has  been  entered.  Nor  can  the  user  re-open  (and  modify)  existing  input  files. 
The  relatively  small  input  data  requirements  did  not  justify  the  complexities  inherent  in 
adding  this  type  of  procedural  flexibility. 

•  While  the  output  does  provide  the  most  relevant  information  in  a  straightforward  manner, 
work  print-outs  are  not  particularly  aesthetic. 

•  Where  the  user  inputs  an  invalid  data  entry  (e.g.  inputting  an  alpha  character  where  a 
numeric  value  is  required),  the  program  will  immediately  terminate  (i.e.,  "bomb").  Again, 
this  was  not  regarded  as  a  significant  deficiency  as  the  limited  amount  of  required  input 
data  generally  makes  re-inputting  information  simple  and  timely. 

•  Where  terminating  an  input  file  is  desirable,  the  user  can  type  the  DOS  command 

contl  -  c 

•  Where  the  user  responds  to  certain  required  input  data  (e.g.  work  zone  length)  by  hitting 
<retum>  without  actually  entering  information,  the  prompt  will  persist  until  data  is  input. 
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B.7  Work  Source  Code 


For  the  purpose  of  "historical  record",  the  work  source  code  is  presented  on  the  following 


pages. 


DIMENSION  QI(30) ,QO(30) ,DI(30) ,DO(30) ,AI(30) ,AO(30) 

real      egi,ego 

real    ila, ild, il , i2a, i2d, i2 , vl, v2 

REAL   TDI (30), TDO (30) 

real   mu,wlength, lsistd, lsostd, vd, lo, li 

character    *10,    name 

character    *10,    street (2) 

character    *  1 , opt 

character    *l,optl 

character  *30 , jobtype, comment 

integer  job,sttype 

REAL  LSI,LSO,lgi,lgo 

print  *, 'name  of  highway  (up  to  10  characters) ' 

read  8 9 7, name 
897      format  (alO) 
71       format (a30) 

print  *,'if  the  highway  is  east-west,  type  1' 

print  *, 'if  the  highway  is  north-south,  type  2' 

read  * , sttype 

if (sttype.eq. 1)  then 

street (1 ) = ' eastbound' 

street (2 ) = 'westbound' 

endif 

if (sttype.eq. 2)  then 

street ( 1 )  =  ' northbound ' 

street (2 ) = ' southbound' 

endif 

print  *, 'hours  of  analysis' 

read  * , n 

print  *, 'maximum  green  time  for  ', street (1) , '  direction  (mins) ' 

read  * ,  gmi 

print  *, 'maximum  green  time  for  ', street (2 ), '  direction  (mins)' 

read  * ,  gmo 

gmi=gmi*60 . 

gmo=gmo*60 . 

print  *, 'work  zone  length  (miles) ' 

read  *,wlength 

wlength=wlength*52  80 . 0 

t=3600. 

if (sttype. eq.l)  then 

print  *, 'vehicle  demand  rate-eastbound  (veh/hr) ' 

do  100  i=l,n 

print  * , 'hour' , i 

read  *,qi(i) 
100      continue 

print  *, 'vehicle  demand  rate-westbound  (veh/hr)' 

do  101  i=l,n 
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68  print 
print 
print 
print 


print  *,'hour',i 
read  * , qo ( i ) 
101      continue 
endif 

if (sttype.eg. 2 ) then 

print  *, 'vehicle  demand  rate-northbound  (veh/hr) ' 
do  890  i=l,n 
print  *,'hour',i 
read  * , gi ( i ) 

890  continue 

print  *, 'vehicle  demand  rate-southbound  (veh/hr)' 
do  891  i=l,n 
print  * , 'hour' , i 
read  *,qo(i) 

891  continue 
endif 

print  *,  'speed  limit   or   estimated  prevailing  speed' 
print  * , ' in  the  absence  of  a  work  zone  (mile/hr) ' 

read  * ,  vd 

vd=vd*5280./3600. 

print  *,  'job  type' 

'Asphalt  Concrete  Overlay      =  1' 
'Chip  Seal  =  2' 

'other  if  pilot  car  is  not  used  =  3' 
'other  if  pilot  car  is  used     =  4' 

read  *,  job 

if (job. eg. 1)  jobtype= 'Asphalt  Concrete  Overlay' 

if (job. eg. 2)  jobtype= 'Chip  Seal' 

if(job.eg.3)  jobtype= 'Other   without  Pilot  Car' 

if(job.eg.4)  jobtype= 'Other   with  Pilot  Car' 

67       print  *, 'Comment  (up  to  30  characters)' 

read  71, comment 

write(2,879) 
879      format (///  '  ANALYSIS  RESULTS',/, 

',//) 

write (2 , 72 )  name, n, jobtype, wlength/5280 .  ,  street (1 ) , gmi/60 . 

*  street (2) ,gmo/60. , vd*3  600. /5280. , comment 

72     format (6x, 'Highway' , 31x,al0, //, 6x, 'Hours  of  Analysis ', 18x, i2 , 

*  4x, 'hour (s) ',//, 6x, 'Job  Type' , 17x, a30 , //6x, 'Work  Zone  Length' 
*,  15x, f6.2, 4x,  'mile(s)  ' , //, 6x,  'Maximum  Green  Time' , /, 10x,al0, 
*18x, f7 .2, 4x, 'min' , / , lOx, alO, 18x, f 7 . 2 , 4x, 'min' , //, 6x, 

* 'Normal  Speed' ,21x, f6.2,4x, 'mph' , //, 6x 
* 'Comments    :',10x,a30) 

write(2,374)  street (1) , street (2 ) 
374  format(////, 

*'  Vehicle  Demand  Rate  (veh/hr) ',//, llx, 9x, 

*  al0,5x,al0) 

do  32  kk=l,n 

write (2, 351)  kk, gi (kk) , go (kk) 
351       format (6x, 'hour' , i2,2f 15.0) 
32       continue 

criteria=12 . 21 
if (job. eg. 1)  goto  701 
if (job. eq. 2)  goto  702 
if(job.eq.3)  goto  703 
if(job.eg.4)  goto  704 


66 


701       ssi=1083.803 
sso=1083.803 
st=1083.803 
ila=1.0 
ild=23.19 
il=ila+ild 
i2a=1.0 
i2d=23.19 
i2=i2a+i2d 
vl=0.0 
v2=0.0 
vi  =  33. 74384 
vo=33.  74384 
spdi=33. 74384 
spdo=33. 74384 
spd=33. 74384 
lsi=30.1473 
lso=30.1473 
sistd=188.5843 
sostd=188.5843 
vistd=10. 49788 
vostd=10. 49788 
lsistd=0.0 
lsostd=0.0 
lgi=13.5503 
lgo=13.5503 
spdivar=vistd**2 
spdovar=vostd**2 
ski=ssi/3600 
sko=sso/3600 
li=30.1473 
lo=30.1473 
goto  800 


702       ssi=1018.4893 
sso=1018.4893 
st=1018.4893 
ila=1.0 
ild=3.89 
il=ila+ild 
i2a=1.0 
id2=3.89 
i2=i2a+i2d 
vl=0.0 
v2=0.0 

vi=29. 076315 
vo=29. 076315 
spdi=29. 076315 
spdo=29. 076315 
spd=29. 076315 
lsi=23.32 
lso=23.32 
sistd=302. 93275 
sostd=302. 93275 
vistd=7. 67177 
vostd=7. 67177 
lsistd=0.0 
lsostd=0.0 
lgi=6.9618 
lgo=6.9618 
spdivar=vistd**i 
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spdovar=vostd**2 

ski=ssi/3600 

sko=sso/3600 

li=23.32 

lo=23.32 

goto  800 

703       ssi=1090.58 
sso=1090.58 
st=1090.58 
ila=1.0 
ild=15.91 
il=ila+ild 
i2a=1.0 
i2d=15.91 
i2=i2a+i2d 
vl=0.0 
v2=0.0 
vi=44.76 
vo=44 .76 
spdi=44 .76 
spdo=44.76 
spd=44.76 
lsi=34. 97136 
lso=34. 97136 
sistd=380.6679 
sostd=380. 6679 
vistd=10. 97534 
vostd=10. 97534 
lsistd=104.9463 
lsostd=104.9463 
lgi=6.9618 
lgo=6.9618 
spdivar=vistd**2 
spdovar=vostd**2 
ski=ssi/3600 
sko=sso/3600 
li=34. 97136 
lo=34. 97136 
goto  800 


704      ssi=1047.679 
sso=1047.679 
st=1047.679 
ila=1.0 
ild=2.10 
il=ila+ild 
i2a=1.0 
i2d=2.10 
i2=i2a+i2d 
vl=0.0 
v2=0.0 
vi=29.75 
vo=29.75 
spdi=29.75 
spdo=29.7  5 
spd=29.75 
lsi=15.37 
lso=15.37 
sistd=160.1344 
sostd=160.1344 
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vistd=3. 233655 

vostd=3. 233655 

lsistd=0. 0 

lsostd=0 . 0 

lgi=8.22 

lgo=8.22 

spdivar=vistd**2 

spdovar=vostd*  *2 

ski=ssi/3600 

sko=sso/3600 

li=15.37 

lo=15.37 

goto  800 

800      DO  16  1=1, N 

QI(I)=QI(I) /3600. 

QO(I)=QO(I)/3600. 
16       continue 

Tdelay=wlength/spd-wlength/vd 
St=St/3600. 

C=GMI+GMO+2 . * WLENGTH/ SPD+LSO+LSI 
SI=St*GMI/C 
SO=St*GMO/C 
1  =  0 
1        IF(I.GE.N)  GOTO  9999 
AI(I+1)=0.0 
AO(I+1)=0.0 
1=1+1 

GI  =  01 (I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgo)-QO(I)* (lgo-lgi) ) / 

*  (St*(St-(QI(I)+QO(I))))+lgi 

GO  =  QO(I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgi)-QI(I)*(lgi-lgo) ) / 

*  (St*(St-(QI(I)+QO(I) ) ) )+lgo 

IF( (GI.LE.GMI.AND.GI.GE.0.0) .AND. (GO . LE. GMO. AND. GO .GE . 0 . 0 ) ) THEN 

c        calculating   expected  travel  time 

888        timei  =  wlength* (spdi**2+spdivar) /spdi**3 
timeo  =  wlength* (spdo**2+spdovar ) /spdo**3 

c        calculating  variance  of  travle  time 

timeivar=wlength**2* ( (spdi**2+3*spdivar) /spdi**4 

*  - ( (spdi**2+spdivar) /spdi**3) **2) 
timeovar=wlength**2* ( (spdo**2+3*spdovar) /spdo**4 

*  -( (spdo**2+spdovar) /spdo**3)**2) 

c        P(headway>C) 

pci=exp (-1. *qi (i) * criteria) 
pco=exp ( -1 . *qo (i) * criteria) 

c        expected  extended  green  time 

egi= (1-exp (-l*qi (i) * criteria) ) / (exp(-l*qi (i) * criteria) ) /qi (i) 
ego= (1-exp ( -l*qo ( i) * criteria) ) / (exp ( -l*qo (i) * criteria) ) /qo (i) 

c        variance  of  extended  green  time 
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egivar=(l-exp(-l*qi (i) *criteria) ) / (exp(-l*qi (i) * 

criteria) ) /qi (i) **2+ (l-exp(-l*qi (i) *criteria) ) 
/ (exp(-2*qi(i) *criteria) ) /qi(i)**2 

egovar= (l-exp(-l*qo (i) * criteria) ) / (exp ( -l*qo (i) *criteria) ) 
/qo(i) **2+(l-exp(-l*qo(i) *criteria) ) / (exp ( -2*qo (i) * 
criteria) ) /qo(i) **2 

expected  green  time 

gi=qi (i) /ski/ (1-qi (i) /ski-qo (i) /sko) * (timei+  timeo+qo (i) / 

sko*egi+ (1-qo (i) /sko) *ego+li+lo) 
go=qo(i) /sko/ (l-qo(i) /sko-qi (i) /ski) * (timeo+  timei+qi (i) / 

ski*ego+ (1-qi (i) /ski) *egi+lo+li) 


variance  of  green  time 

xo=i2*qo (i) * (timeo+gi+egi+timei+li+lo) / 

( sko* *2* (1-qo (i) /sko) **3) 
yo=v2/ (sko/qo (i ) ) **2 
zo=qo(i)**2/ ( (sko-qo(i) ) **2) 
xi=il*qi (i) * (timei+go+ego+timeo+lo+li ) / 

(ski**2* (1-qi (i) /ski) **3) 
yi=vl/(ski/qi(i) )**2 
zi=qi(i)**2/ (ski-qi (i) )**2 

aoo=xo+yo+zo* (timeovar* (ski+qi (i) ) / (ski-qi (i) ) +xi+yi+zi* 
( timeivar * ( sko+qo ( i ) ) 

/ (sko-qo (i ) ) +egovar+timeovar) +egivar+timeivar) 
govar=aoo/ (l-zi*zo) 

aoi=xi+yi+zi* (timeivar* (sko+qo (i) ) / (sko-qo (i) ) +xo+yo+zo 
* ( timeovar* ( ski+qi ( i ) ) 

/ (ski-qi (i) ) +egivar+timeivar) +egovar+timeovar) 

givar=aoi/ (l-zo*zi) 

expected  red  time 

ri=go+timei+ego+timeo+li+lo 
ro=gi+timeo+egi+timei+lo+li 

variance  of  red  time 

rivar= (sko+qo (i) ) / (sko-qo (i) ) *timeivar+govar+egovar+timeovar 
rovar= (ski+qi (i) ) / (ski-qi (i) ) *timeovar+givar+egivar+ timeivar 

total  delay  per  cycle 

totali=qi (i) / (2* (1-qi (i) /ski) ) * (ri**2+rivar+il*ri/ 
(ski* (1-qi (i) /ski) ) +vl/ (ski*qi (i) ) ) 

totalo=qo(i) / (2* (l-qo(i) /sko) ) * (ro**2+rovar+i2*ro/ 
(sko* (1-qo (i) /sko) ) +v2/ (sko*qo(i) ) ) 

average  delay 

avgi=totali/ ( (timei+li+go+ego+timeo+gi+egi+lo) *qi (i) ) 
avgo=totalo/ ( ( timei+li+go+ego+timeo+gi+egi+lo) *qo (i) ) 
di ( i ) =avgi+tdelay 
do ( i ) =avgo+tdelay 
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total  delay  for  the  hour 


199 
899 

799 


757 


tdi ( i ) =di ( i ; 
tdo(i)=do(i! 


rt*qi(i) 
rt*qo(i) 


write(2,199)  i 

format (///// '  hour',i2,'   under-saturated') 

format (///// '  hour',i2,'   over-saturated') 

write (2 ,799)  street (1) ,ai(i) , ai(i+l) , street (2) , ao ( i) , ao ( i+1) 
format (//'  Initial  Queue      Residual  Queue',//, 

2x,al0,7x, f6.0,10x,f6.0, / , 2x, alO, 7x, f 6 . 0, lOx, f6.0) 
write (2, 757)  street (1) , di(i) /60. ,tdi(i) /3600.  , 
street (2) ,do(i) /60. ,tdo(i) /3  6  00. 

f ormat (//, 5x, 'Delay   -   Average  (min/veh)     Total  (veh-hr) ' 
, //,2x,al0,7x, f 6. 1,1 Ox, f 7 . 2 , /,2x,al0,7x, f 6 . 1 , lOx,  f 7 . 2 ) 
GOTO  1 
ENDIF 


599 


299 


257 


IF( (GI.GT.GMI.OR.GI.lt. 0.0) -AND. (GO.GT .GMO .OR. GO . LT. 0.0) )  THEN 

C=GMI+GMO+2  *WLENGTH/SPD+LSO+LSI 

SI=St*GMI/C 

SO=St*GMO/C 

AI(I+1)=QI(I) *T-SI*T 

AO(I+l)=QO(I)*T-SO*T 


IF(AI(I+1) .LE.0.0 
IF(AI(I+1) .GT.0.0 
IF(AI(I+1) .LE.0.0 


.AND. 

AO(I+l) 

.  LE 

.AND. 

AO(I+l) 

.LE 

.AND. 

AO(I+l) 

.GT 

0.0)  GOTO  21 
0.0)  GOTO  22 
0.0)  GOTO  23 


TDI (I) =AI (1+1) *T*0 . 5+tdelay*si*t+st* (c-gmi] 
TDO(I)=AO(I+l) *T*0.5+tdelay*so*t+st* (c-gmo) 
di(i)=tdi(i) / (ai (i) +qi (i) *t) 
do(i)=tdo(i) / (ao(i)+qo(i) *t) 
write(2,599)  i 
f ormat (///// ' 


tgmi*3600/c 
tgmi*3600/c 


hour',i2,'   over-saturated') 

street (1) ,ai(i) , ai (i+1) , street (2)  ,  ao (i) , ao (i  +  1) 

Initial  Queue      Residual  Queue' , // 
0,10x,f6.0, /,2x,al0,7x,f6.0,10x, f6.0) 
street (1) ,di (i) /60 . , tdi (i) /3600 . , 
street (2) ,do(i) /60. , tdo(i) /3  600. 
format (//, 5x, 'Delay   -   Average  (min/veh) 
//,2x,al0,7x,f6.1,10x,f7.2, / , 2x, alO , 7x, f 6 . 1 
GOTO  2 


write(2,299) 
format (//' 
2x,al0,7x, f6 
write (2, 257 


Total  (veh-hr) 
,10x,f7.2) 


IF(I.GE.N) 
1  =  1  +  1 


GOTO  9999 


AI ( 1+1 ) =AI ( I ) +QI ( I ) *T-SI *T 
AO ( 1+1 ) =AO ( I ) +QO ( I ) *T-SO*T 

IF(AI(I+1) .GT.0.0.AND.AO(I+l) .GT.0.0)  THEN 

TDI(I)=(AI(I)+AI(I+1) )*T*0.5+tdelay* (si*t) 
+st* (c-gmi) *0 . 5*gmi*3  600/c 

TDO(I)=(AO(I)+AO(I+l) ) *T*0.5+tdelay*(so*t) 
+st* (c-gmo) *0 . 5*gmo*3  600/c 

di ( i ) =tdi ( i ) / (ai ( i ) +qi ( i ) *t ) 

do(i)=tdo(i) / (ao(i)+qo(i)*t) 

write(2,899)  i 
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write (2, 7 99)  street (1) ,ai(i) ,ai(i+l) , street (2) ,ao(i) ,ao(i+l) 
write (2, 757)  street (1) ,di(i)/60. , tdi(i) /3  600.  , 
*  street (2) , do(i) /60. , tdo(i) /3  600. 

GOTO  2 

ENDIF 

21       IF(AI (1+1) .LE.0.0.AND.AO(I+l) . LE . 0 . 0 )  THEN 

ai(i+l)=0.0 
ao(i+l) =0. 0 
TI=AI(I)/(SI-QI(I) ) 
TO=AO ( I ) / ( SO-QO ( I ) ) 

IF(TI-TO)  3,  4,  5 
3  GI1=(QI(I) * (2*WLENGTH/SPD+GMO+LSI+LSO) ) /(St-QI(I) ) +lgi 

Cl=GIl+GMO+2*WLENGTH/SPD+LSO+LSI 
C=GMI+GMO+2*WLENGTH/SPD+LSO+LSI 
S01=GMO*St/Cl 

T02=(AO(I) -SO*TI+S01*TI) / (SOl-QO(I) ) 
GI2=QI (I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgo) -QO (I ) * (lgo-lgi) ) 

/(St*(St-(QI(I)+QO(I))))+lgi 
G02=QO(I) * (St*(2*WLENGTH/SPD+LSO+LSI+lgi)-QI(I)*(lgi-lgo) ) 

/(St*(St-(QO(I)+QI(I) ) ) )+lgo 
C2=GI2+G02+2*WLENGTH/SPD+LSO+LSI 

TDI(I)=AI(I)  *TI*0.5  +  st*(c-grni)  *0  .  5*gmi*ti/c 

*  +0.5*(l-GIl/Cl)**2*Cl/(l-QI(I)/St)*QI(I)*(TO2-TI) 

*  +0.5* (1-GI2/C2) **2*C2/ (1-QI (I) /St) *QI (I) * (T-T02) 

*  +tdelay* (si*ti+gi(i) * ( to2-ti) +qi (i) * (t-to2) ) 

TDO(I)=0.5* (AO(I)+AO(I)+QO(I)*TI-SO*TI) *TI 

*  +0.5*(AO(I)+QO(I) *TI-SO*TI)*(T02-TI) 
+0.5*(l-GO2/C2)**2*C2/(l-QO(I)/St)*QO(I)*(T-TO2) 
+tdelay* (so*ti+sol* (to2-ti) +qo (i) * (t-to2) ) 

*  +st* (c-gmo) *0 . 5*gmo*ti/c 

*  +st* (cl-gmo) *0.5*gmo* (to2-ti) /cl 
di(i)=tdi(i)/(ai(i)+qi(i)*t) 
do(i)=tdo(i) / (ao(i)+go(i)*t) 

write(2,899)  i 

write (2, 7 99)  street (1) , ai (i) , ai (i+1) , street (2) , ao (i) , ao (i+1) 

write (2, 757)  street (1) ,di(i)/60. , tdi(i) /3600.  , 

*  street (2) ,do(i) /60. , tdo(i) /3600. 

GOTO  1 

G01=(QO(I) * (2*WLENGTH/SPD+GMI+LSO+LSI) ) / (St-QO(I) ) +lgo 

Cl=G01+GMI+2*WLENGTH/SPD+LSI+LSO 

C=GMO+GMI+2*WLENGTH/SPD+LSI+LSO 

SIl=GMI*St/Cl 

TI2=(AI(I)-SI*TO+SIl*TO)/(SIl-QI(I) ) 

G02=QO(I) * (St* (2*WLENGTH/SPD+LSO+LSI+lgi) -QI (I) * (lgi-lgo) ) 

*  /(St*(St-(QI(I)+QO(I) ) ) )+lgo 
GI2=QI (I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgo) -QO(I) * (lgo-lgi) ) 

*  /(St*(St-(QO(I)+QI(I) ) ) )+lgi 
C2=GI2+G02+2*WLENGTH/SPD+LSO+LSI 


* 


TDO(I) =AO(I) *TO*0 . 5+st* (c-gmo) *0 . 5*gmo*to/c 

M*QO(I)*(TI2-TO) 
)(I) * (T-TI2) 


)=AO(I) *TO*0.5+st* (c-gmo) *0 . 5*gmo*to/c 
+0.5* (1-GOl/Cl) **2*Cl/(l-QO(I) /St)*QOI 
+0.5* (1-G02/C2) **2*C2/(l-QO(I) /St)*QOl 
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+tdelay* (so*to+qo(i) * ( ti2-to) +qo (i) * (t-ti2) ) 

TDI(I)=0.5* (AI(I)+AI(I)+QI(I)*TO-SO*TO) *T0 
+0 . 5* (AI (I ) +QI (I ) *TO-SI*TO) * (TI2-TO) 
+0.5* (1-GI2/C2) **2*C2/(1-QI(I) /St)*QI(I) * (T-TI2) 
+tdelay* (si*to+sil* (ti2-to) +qi (i)*(t-ti2) ) 
+st* (c-gmi) *0 .5*gmi*to/c 
+st* (cl-gmi) *0.5*gmi* (ti2-to) /cl 

di ( i ) =tdi ( i ) / ( ai ( i ) +qi ( i ) * t ) 

do(i)=tdo(i) / (ao(i)+qo(i)*t) 

write(2,899)  i 

write (2,799)  street (1) , ai (i) , ai (i  +  1) , street (2)  ,  ao (i) , ao (i  +  1) 

write (2, 757)  street (1) ,di(i) /60. ,tdi(i) /3  600.  , 
street (2) ,do(i) /60. , tdo(i) /3  600. 

GOTO  1 


C=GMI+GMO+2  *WLENGTH/SPD+LSO+LSI 

GI1=  QI (I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgo) -QO(I) * (lgo-lgi) ) 

*  /(St*(St-QI(I)+QO(I) ) )+lgi 

G01=  QO(I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgi) -QI (I) * (lgi-lgo) ) 

*  / (St* (St-QO(I)+QI(I) ) )+lgo 
Cl=   GIl+G01+LSO+LSI+2*WLENGTH/SPD 

TDI(I)=AI(I) *TI*0.5+st* (c-gmi) *0 . 5*gmi*ti/c 

*  +0.5* (1 -Gil /CD  **2*C1/(1-QI(I) /St) *QI(I) * (T-TI) 

*  +tdelay* (si*ti+qi (i) * (t-ti) ) 
TDO(I)=AO(I) *TO*0.5+st* (c-gmo) *0 . 5*gmo*to/c 

*  +0.5* (1-G01/C1)**2*C1/(1-Q0(I) /St) *QO(I) * (T-TO) 

*  +tdelay* (so*to+qo (i) * ( t-to) ) 

di(i)=tdi(i) /(ai(i)+qi(i)*t) 
do(i)=tdo(i) / (ao(i) +qo(i) *t) 

write(2,899)  i 

write (2 , 799)  street (1) ,ai(i) ,ai(i+l) , street (2) , ao (i) , ao (i+1) 

write (2, 7 57)  street (1) ,di(i) /60. , tdi(i) /3  600. , 

*  street (2) ,do(i) /60. , tdo(i) /3  600. 

GOTO  1 
ENDIF 


22         IF(AI(I+1) .GT.0.0.AND.AO(I+l) .LE.0.0)  THEN 
ao(i+l) =0. 0 
TO=AO(I) /(SO-QO(I) ) 

G01=QO(I) * (2*WLENGTH/SPD+GMI+LSI+LSO) / (St-QO(I) ) +lgo 
Cl=G01+GMI+2*WLENGTH/SPD+LSO+LSI 
C=GMO+GMI+2*WLENGTH/SPD+LSO+LSI 

SIl=GMI*St/Cl 

SI=GMI*St/C 

SO=GMO*St/C 

AI(I+1)=AI(I)+QI(I)*T-SI*T0-SI1* (T-TO) 
IF(AI(I+1) .GT.0.0)  THEN 


73 


TDI(I)=(AI (I)+AI(I)+QI(I) *TO-SI*TO) *0.5*TO 

*  +(AI(I)+QI(I) *T0-SI*T0+AI(I+1) ) *0.5* (T-TO) 

*  +tdelay* (si*to+sil* (t-to) ) 

*  +st* (c-gmi) *0 . 5*gmi*to/c 

*  +st*(cl-gmi) *0.5*gmi* (t-to) /cl 
TDO(I)=AO(I) *TO*0.5 

*  +0.5* (1-G01/C1)**2*C1/ (l-QO(I)/St)*QO(I)*(T-TO) 

*  +tdelay* (so*to+qo(i) * (t-to) ) 

*  +st* (c-gmo) *0 . 5*gmo*to/c 

di(i)=tdi(i) / (ai (i) +gi (i) *t ) 
do(i)=tdo(i) / (ao(i)+go(i) *t) 
AO(I+1)=0.0 

write(2,899)  i 

write (2, 799)  street (1) , ai (i) , ai (i+1) , street (2) , ao (i) , ao (i+1) 

write (2, 757)  street (1) ,di(i) /60. ,tdi(i) /3  600. , 

*  street (2) , do(i) /60. , tdo(i) /3  600. 


GOTO  2 

ENDIF 

IF(AI(I+1) .LE.0.0)  THEN 

ai(i+l)=0.0 
if (ai (i) .eq.O.O.and.ao(i) .eq.0.0)  goto  888 

T01=(AI(I)-SI*TO+SIl*TO) /(SIl-QI(I) ) 
GI2=QI (I) * (St* (2*WLENGTH/SPD+LSI+LS0+lgo) -QO (I ) * ( lgo-lgi ) ) / 
(St*(St-(QI(I)+QO(I) ) ) )+lgi 
G02=QO(I) *(St* (2*WLENGTH/SPD+LSO+LSI+lgi)-QI(I)*(lgi-lgo) ) / 

*  (St*(St-(QO(I)+QI(I) ) ) )+lgo 

C2=GI2+G02+2*WLENGTH/SPD+LS0+LSI 

TDI ( I )  =  ( AI ( I ) +AI ( I ) +TO*QI ( I ) -TO*SI ) *  0 . 5*TO 

*  +(AI(I)+TO*QI(I)-SI*TO)*(T01-TO) *0.5 
+0.5*(1-GI2/C2)**2*C2/ (1-QI (I) /St) *QI (I) * (T-TOl) 
+tdelay* (si*to+sil* (tol-to) +gi (i) * (t-tol) ) 
+st* (c-gmi ) *0 . 5*gmi*to/c 
+st* (cl-gmi) *0. 5*gmi* (tol-to) /cl 

TDO(I)=AO(I) *TO*0.5+st* (c-gmo) *0 . 5*gmo*to/c 

+0.5* (1-G01/C1) **2*C1/(1-Q0(I) /St)*QO(I) *(T01-TO) 
+  0 . 5  * ( 1 -G02 /C2 ) *  *2  *C2 / ( 1 -QO ( I ) /St ) *QO (I) * ( T-TOl ) 
+tdelay* (so*to+qo(i) * (tol-to) +qo (i) * (t-tol) ) 


* 


di ( i ) =tdi (i)/(ai(i)+qi(i)*t) 

do ( i ) =tdo ( i ) / ( ao ( i ) +qo ( i ) *  t ) 
write(2,899)  i 

write (2, 799)  street (1) , ai (i) , ai (i+1) , street (2) , ao(i) ,ao(i+l) 
write (2, 7 57)  street (1) ,di(i) /60. ,tdi(i) /3600. , 
*  street (2) ,do(i)/60. , tdo(i) /3  600. 

GOTO  1 

ENDIF 
ENDIF 

23         IF(AO(I  +  l)  .GT.0.0.AND.AKI  +  1)  .LE.0.0)  THEN 
ai(i+l)=0.0 

TI=AI(I)/(SI-QI(I) ) 

GI1=QI (I) * (2*WLENGTH/SPD+GMO+LSO+LSI) / (St-QI (I) ) +lgi 

Cl=GIl+GMO+2*WLENGTH/SPD+LSI+LSO 
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C=GMI+GMO+2*WLENGTH/SPD+LS0+LSI 

S01=GM0*St/Cl 

SO=GMO*St/C 

SI=GMI*St/C 

A0(I+1)=A0(I)+Q0(I) *T-S0*TI-S01*<T-TI) 
IF(A0(I+1) .GT.O.O)  THEN 


TDO ( I) =(AO (I) +AO (I) +00 (I) *TI-SO*TI) *0.5*TI 

*  +(AO(I)+QO(I) *TI-SO*TI+AO(I+l) ) *0.5* (T-TI) 

*  +tdelay* (so*ti+sol* (t-ti) ) 

*  +  st* (c-gmo) *0 . 5*gmo*ti/c 

*  +st* (cl-gmo) *0 . 5*gmo* (t-ti) /c 
TDI(I)=AI(I) *TI*0.5+st*(c-gmi) *0 . 5*gmi*ti/c 

*  +0.5* (1-GI1/C1)**2*C1/ (l-QI(I) /St) *QI(I)* (T-TI) 

*  +tdelay* (si*ti+qi (i) * (t-ti) ) 
di(i)=tdi(i) / (ai (i) +qi (i) *t) 
do(i)=tdo(i) / (ao(i)+qo(i) *t) 

AI(I+1)=0.0 

write (2, 899)  i 
write (2, 7 99)  street (1) ,ai(i) , ai (i+1) , street (2) , ao ( i) , ao (i+1) 
write (2 ,757)  street (1) ,di(i) /60. , tdi(i) /3600. , 

*  street (2) ,do(i) /60. , tdo(i) /3  600. 


GOTO  2 
ENDIF 

IF (AO ( 1+1 ) . LE . 0 . 0 )  THEN 
ao(i+l)=0.0 
if  (aid)  .eq.O.O.and.ao(i)  .eq.0.0)  goto  888 

TI1= (AO(I) -SO*TI+S01*TI) / (SOl-QO(I) ) 
G02=QO(I) * (St* (2*WLENGTH/SPD+LSO+LSI+lgi) -QI (I) * (lgi-lgo) ) / 
(St*(St-(QO(I)+QI(I) ) ) )+lgo 

GI2=QI (I) * (St* (2*WLENGTH/SPD+LSI+LSO+lgo) -QO(I) * (lgo-lgi) ) / 
(St*(St-(QI(I)+QO(I) ) ) )+lgi 

C2=G02+GI2+2*WLENGTH/SPD+LSI+LSO 

TDO(I)=(AO(I)+AO(I)+TI*QO(I) -TI*SO) *0.5*TI 
+(AO(I)+TI*QO(I)-SO*TI) * (TI1-TI) *0.5 

*  +0.5* (l-G02/C2)**2*C2/(l-QO(I) /St)*QO(I)* (T-TI1) 

*  +tdelay* (so*ti+sol* (til-ti) +qo (i) * (t-til) ) 

*  +st* (c-gmo) *0 . 5*gmo*ti/c 

*  +st* (cl-gmo) *0 .5*gmo* (til-ti) /c 
TDI(I)=AI(I)*TI*0.5+st*(c-gmi)*0.5*gmi*ti/c 

*  +0.5* (1-GI1/C1)**2*C1/(1-QI(I) /St) *QI(I)* (TI1-TI) 

*  +0.5* (1-GI2/C2) **2*C2/(1-QI(I) /St)*QI(I)* (T-TI1) 

*  +tdelay* (si*ti+qi(I) * ( til-ti )+qi (i) * (t-til) ) 


* 


di(i)=tdi(i) /(ai(i)+qi(i)*t) 

do (i ) = tdo ( i ) / ( ao ( i ) +qo ( i ) *  t ) 
write(2,899)  i 
write  (2,  7  99)  street  (1)  , aid)  ,  ai  (i  +  1)  ,  street  (2)  ,ao(i)  ,ao(i  +  i; 
write (2, 7 57)  street (1) ,di(i)/60. ,tdi(i) /3  600. , 
street (2) ,do(i) /60. , tdo(i) /3  600. 

GOTO  1 
ENDIF 
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ENDIF 
ENDIF 
IF( (GI.GT.GMI.OR.GI.LT. 0.0) .AND. (GO. LE.GMO . AND.GO.GE. 0 . 0 ) )  THEN 

GO=QO(I) * (2*WLENGTH/SPD+GMI+LSI+LSO) / (St-QO(I) ) +lgo 

Cl=  GMI+GO+2*WLENGTH/SPD+LSO+LSI 

AI(I+1)=QI(I) *T-St*GMI/Cl*T 

AO(I+1)=0.0 

GI=GMI 

IF  (AI(I+1) .LE.0.0. AND. AO(I+l) .LE.0.0)  GOTO  888 
IF  (AI(I+1) .GT.0.0.AND.AO(I+l) -LE.0.0)  GOTO  22 
DI(I)=0.5*AI(I+1) *T/ (QI(I)*T) 

DO(I)=0.5*(1-GO/C1) **2*Cl/(l-QO(I) /St) 

write(2,899)  i 
write (2, 799)  street (1) , ai (i ) , ai (i+1) , street (2) , ao (i) , ao (i+1 ) 
write (2, 757)  street (1) ,di(i)/60. , tdi (i) /3600 . , 

*  street (2) ,do(i) /60. , tdo(i) /3  600. 

GOTO  2 

ENDIF 

IF ( ( GO . GT . GMO . OR . GO . LT . 0 . 0 ) . AND . ( GI . LE . GMI . AND . GI . GE . 0 . 0 ) )  THEN 

GI=QI ( I ) * ( 2  *WLENGTH/SPD+GMO+LSO+LSI ) / ( St-QI ( I ) ) +lgi 

Cl=  GMO+GI+2*WLENGTH/SPD+LSI+LSO 

AO(I+l)=QO(I) *T-St*GMO/Cl*T 

AI(I+1)=0.0 

GO=GMO 

IF(AO(I+l) .LE.0.0  .AND.  AI ( 1+1 ) . LE. 0 . 0 )  GOTO  888 
IF(AO(I+l) .GT.0.0  .AND.  AI ( 1+1 ) . LE. 0 . 0 )  GOTO  23 
DO(I)=0.5*AO(I+1) *T/ (QO(I)*T) 
DI(I) =0.5* (1-GI/C1)**2*C1/(1-QI(I) /St) 
write (2, 899)  i 
write (2 ,799)  street (1) ,ai(i) ,ai(i+l) , street (2) ,ao(i) ,ao(i+l) 
write (2, 757)  street (1) ,di(i)/60. , tdi (i ) /3600 .  , 

*  street (2) ,do(i)/60. , tdo (i) /3  600 . 

GOTO  2 

ENDIF 
9999  write(2,174)  street ( 1) , street (2 ) 

103  format (///, 

*'     Total  Delay  (veh-hrs) ' , // , 4x, 

*  al0,2x,al0) 
803             format (//, 

*'     Average  Delay  (min/veh) ' , // , 4x, 

*  al0,2x,al0) 

174  format (///////, 

*'  Total  Delay  (veh-hrs) ',//, 4x, 9x, 

*  al0,5x,al0) 
DO  49  1=1, N 
gi(i)=qi(i) *3600. 
qo(i)=qo(i) *3600. 

WRITE (2, 51)  i,TDI(I) /3  600. , TDO ( I ) /3  600. 
51  FORMAT (' hour '  ,i2,2F15.2) 
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49  CONTINUE 

write (2, 69) 
69  format (/////) 

st=st*3600. 

print  *, 'Do  you  want  to  change  options  ?  (Yes  =  y  ,  No  =  n) ' 
read  61,  opt 
61  format (al) 

i  f ( opt . ne . ' y ' . and . opt . ne . ' n ' )  then 

75  print  74,  opt 

74  format (2x, al, 4x, ' incorrect  response   ,   Please  try  again') 

read  61, opt 

if (opt .ne. 'y' .and. opt .ne. 'n' )  goto  75 
endif 

if  (opt.eq.'n')  goto  62 
i=0 

print  *, 'Change  Work  Zone  Length  ?  (Yes  =  y,  No  =  n) ' 
read  61,  optl 
if (optl .ne . 'y' . and. optl .ne. 'n' )  then 

76  print  74,  optl 
read  61, optl 

if (optl .ne . 'y' .and. optl .ne . 'n' )  goto  7  6 

endif 

if  (optl.eq. 'y' )  then 

print  *, 'New  Work  Zone  Length  (miles) ' 

read  * ,  wlength 

wlength=wlength*5280 . 

endif 

print  *, 'Change  Maximum  Green  Time  ?  (Yes  =  y,  No  =  n) ' 

read  61, optl 

if (optl .ne. 'y' .and. optl .ne. 'n' )  then 

77  print  74,  optl 
read  61, optl 

if (optl .ne. 'y' . and. optl .ne. 'n' )  goto  77 
endif 

if (optl.eq. 'y' )  then 
print  *, 'New  maximum  green  time  for  ', 

*  street (1),'  direction  (mins)' 
read  * , gmi 

print  *, 'New  maximum  green  time  for  ', 

*  street (2),'  direction  (mins)' 
read  * , gmo 

gmi=gmi*  60 . 

gmo=gmo*60 . 

endif 
print  *, 'Change  Speed  Limit  or  Estimated  Prevailing  Speed' 
print  * , ' in  the  Absence  of  a  Work  Zone  ?  (Yes  =  y,  No  =  n) ' 
read  61,  optl 

if (optl .ne. 'y ' . and. optl .ne . 'n' )  then 

78  print  74,  optl 
read  61, optl 

if (optl .ne . 'y' .and. optl .ne. 'n' )  goto  7  8 

endif 
if (optl.eq. 'y' )  then 

print  * , ' New  Speed  Limit  or  Estimated  Prevailing  Speed' 
print  * , ' in  the  Absence  of  a  Work  Zone  (mile/hr) ' 
read  * ,  vd 
vd=vd*5280./3600. 
endif 

print  *, 'Change  Job  Type  ?  (Yes  =  y,  No  =  n) ' 
read  61, optl 

if (optl .ne. 'y ' .and. optl .ne . 'n' )  then 

77 


79  print  74,  optl 

read  61, optl 

if (optl .ne. 'y' .and. optl .ne . 'n' )  goto  79 
endif 
if  (optl .eq. 'y' )  goto  68 
goto  67 
62  STOP 

END 
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